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I.   INTRODUCTION 

A.   TAYLOR-GORTLER  VORTICES 

Since  the  early  part  of  the  twentieth  century, 
considerable  research  has  shown  that  the  fully  developed 
laminar  flow  along  a  concave  wall  does  not  remain  two- 
dimensional  [Refs.  lf  2,  3],   The  flow  instead  forms  a 
system  of  spiral  vortices,  of  counter  rotating  pairs,  whose 
axes  are  aligned  in  the  principle  flow  direction.   This 
phenomenon  is  the  result  of  the  variations  in  the 
centrifugal  forces  acting  on  the  fluid  particles,  and  is 
known  as  Taylor-Gortler  vortices.   Figure  1  illustrates  the 
type  of  fluid  motion  just  described. 

In  a  channel  that  is  curved  in  the  streamwise  direction, 
those  fluid  particles  located  in  the  center  of  the  flow 
cross-section  are  subject  to  higher  centrifugal  forces  than 
those  fluid  particles  traveling  along  the  channel's  outer 
boundary  wall.   As  a  result,  the  fluid  in  the  center  of  the 
channel  moves  outwardly  toward  the  concave  boundary.   As 
the  process  continues,  the  fluid  particles  near  the  boundary 
wall,  move  in  a  spanwise  direction,  and  finally  radially 
inward  replacing  the  outwardly  moving  particles.   These 
particles  then  come  under  the  same  centrifugal  force  and  the 
process  repeats  itself.   The  resulting  cyclic  motion  causes 
the  formation  of  the  counter  rotating  Taylor-Gortler 
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Primary    Flow,  Direction 


Figure  1.   Schematic  of  Taylor-Gortler  Vortices  in  a  Curved 
Channel . 
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vortices,  considered  to  be  primarily  a  laminar  flow 
phenomenon. 

It  has  been  observed  [Ref.  4]  that  the  heat  transfer 
rate  from  flow  along  a  concave  curved  wall  is  greater  than 
that  for  flow  along  a  straight  wall,  a  phenomenon  attributed 
to  the  additional  mixing  provided  by  the  secondary  motion  of 
the  Taylor-Gortler  vortices. 

There  are  many  possible  applications  that  could  result 
from  a  more  thorough  understanding  of  Taylor-Gortler 
vortices  and  their  effect  on  heat  transfer  and  fluid  flow 
characteristics.   Two  such  applications  could  lead  to 
improved  heat  exchanger  designs  and  improved  turbine  blade 
cooling  [Refs.  5,  6,  7]. 

B.   HISTORY 

The  stability  of  an  inviscid  fluid  flowing  past  a  curved 
boundary,  was  first  considered  by  Lord  Rayleigh  in  1916 
[Ref.  8].   By  assuming  that  the  fluid  was  non-viscous,  he 
determined  that  for  the  motion  to  remain  stable,  its 
circulation  must  increase  with  increasing  radius.   G.  I. 
Taylor,  in  1923,  [Refs.  1,  9],  continued  this  study  with  an 
extensive  analytical  and  experimental  study  of  viscous 
fluids.   His  investigations  focused  on  the  flow  between  two 
cylinders,  in  which  the  inner  cylinder  rotated  while  the 
outer  cylinder  remained  stationary.   Taylor  ascertained  that 
such  flows  become  unstable  when  the  value  of  a  dimensionless 
parameter  exceeded  a  critical  value  of  41.3. 
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The  parameter,  known  as  the  Taylor  number  is  defined  as 


Ta  =  Re  /d/Ri 


where  'd'  is  the  width  of  the  gap,  assumed  small  when 
compared  to  'Ri',  the  radius  of  the  inner  cylinder,  and  'Re' 
is  the  Reynolds  number  based  on  the  peripheral  velocity  of 
the  inner  cylinder.   Taylor  determined  that  for  those  cases 
in  which  the  value  of  the  Taylor  number  exceeded  the 
critical  value,  a  secondary  motion  developed  and  the  Taylor 
vortices  formed.  Figure  2  illustrates  this  fluid  motion. 

Instability  of  a  similar  nature  is  also  observed  when  a 
viscous  fluid  flows  in  a  curved  channel  due  to  a  pressure 
gradient  acting  along  the  channel  wall.   This  problem  was 
first  considered  analytically  by  W.  R.  Dean  [Ref.  10]  in 
1928,  for  a  channel  formed  by  two  concentric  cylinders, 
where  the  radius  of  the  inner  cylinder  was  large  in 
comparison  to  the  spacing  between  the  inner  and  outer 
cylinder  walls.   Dean  concluded  that  there  would  be  an 
initiation  of  flow  instability,  and  the  formation  of 
vortices,  when  a  dimensionless  parameter,  the  Dean  number, 
exceeded  a  value  of  36.   The  Dean  number  is  defined  as: 


De   =   Re 


Ri 


where  'd*  represents  the  channel  half-width,  'Ri •  is  the 
inner  cylinder  radius,  and  'Re*  is  the  Reynolds  number  based 
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Figure  2.   Schematic  of  Taylor  Vortices  between  Cylinders 
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on  the  mean  velocity  of  the  undisturbed  flow  and  the  channel 
half-width  'd'.   The  analytical  work  of  Dean  was  later 
verified  by  W.  H.  Reid  [Ref.  11]/  using  an  approximate 
solution. 

In  1940,  H.  Gortler  [Ref.  2]  studied  the  stability  of 
laminar  boundary  layer  profiles  on  curved  walls  under  the 
influence  of  small  disturbances.   He  found  that  these 
disturbances  were  similar  to  the  vortices  studied  by  G.  I. 
Taylor.   Using  approximate  numerical  calculations,  Gortler 
concluded  that  the  disturbances  of  vortices  were  produced 
only  on  the  concave  boundary  walls  and  that  the  overall  flow 
profile  appeared  to  remain  laminar  in  nature.   These  results 
were  verified,  with  an  exact  solution,  by  G.  Hammerlin,  and 
reported  by  H.  Schlichting  [Ref.  12]  in  1955.   A.  M.  0. 
Smith  completed  an  even  more  extensive  numerical  analysis 
that  further  substantiated  these  findings  [Ref.  3].   The 
results  of  these  numerical  solutions  have  recently  been 
demonstrated  with  the  use  of  hot  wire  anemometry,  laser 
doppler  systems,  and  flow  visualizations  techniques,  [Refs. 
13,  14];  and  in  1976,  Y.  Aihara  [Ref.  15]  conducted  a  non- 
linear analysis  of  Gortler  Vortices. 

As  interest  began  to  develop  concerning  the  effects  of 
the  secondary  flows  associated  with  the  Taylor-Gortler 
vortices,  many  studies  were  published  concerning  the 
influences  of  these  vortices  on  heat  transfer  in  the  laminar 
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and  turbulent  flow  regimes.   In  1955,  F.  Kreith  [Ref.  16], 
studied  the  influence  of  heat  transfer  with  respect  to  the 
curvature  of  the  boundary  wall  for  fully  turbulent  flows. 
He  concluded  that  the  heat  transfer  from  the  heated  concave 
boundary  wall  was  considerably  higher  than  that  transferred 
from  the  convex  boundary  wall  of  the  same  curvature  and 
under  similar  turbulent  conditions. 

In  1965,  L.  Persen  [Ref.  17],  considered  the  special 
cases  of  very  high  and  very  low  Prandtl  number  fluids  and 
related  the  increase  in  heat  transfer  from  a  curved  wall  to 
the  presence  of  the  Taylor-Gortler  vortices. 

There  has  been  only  a  limited  amount  of  published 
literature  dealing  with  the  flow  and  heat  transfer  in  curved 
channels  of  rectangular  cross-section  and  large  aspect 
ratios,  where  the  aspect  ratio  is  equal  to  the  spanwise 
distance  of  the  channel  divided  by  the  channel  height.   Much 
of  what  has  been  published  involves  the  development  of 
numerical  approximations  and  solutions  for  heat  and  mass 
transfer  in  curved  ducts  of  various  geometries.   K.  Cheng 
and  M.  Akiyama  [Ref.  18]  developed  a  numerical  solution  for 
forced  convection  heat  transfer  with  laminar  flows  in  curved 
channels  of  rectangular  cross-section,  but  only  for  small 
aspect  ratios. 

In  1976,  A.  A.  Shibani  and  M.  N.  Ozisik  [Ref.  19],  using 
matched  asymptotic  expansion  techniques  for  a  wide  range  of 
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Prandtl  numbers,  solved  the  heat  transfer  problem  between 
straight  parallel  plates  with  turbulent  flow,  for  the  case 
of  uniform  wall  temperature. 

Other  experimental  and  analytical  studies  that  are  worth 
citing  with  regard  to  this  present  study  are  Y.  Mori  [Ref. 
20],  who  obtained  results  for  hydrodynamically  fully 
developed  flows  with  constant  wall  heat  flux  in  curved 
channels  of  square  cross-section.   W.  M.  Kays  and  E.  Y. 
Leung  [Ref.  21]  reported  solutions  for  turbulent  flow  heat 
transfer  in  a  concentric  circular  tube  annulus  with  a  fully 
developed  velocity  profile  and  constant  heat  rate  per  unit 
length,  for  a  fluid  of  Prandtl  number  0.7.   P.  F.  Brinich 
and  R.  W.  Graham  [Ref.  22]  reported  results  for  turbulent 
flows  in  a  rectangular  curved  channel  with  an  aspect  ratio 
of  6,  for  the  inner  wall  heated,  the  outer  wall  heated,  and 
both  walls  heated. 

In  1973,  at  the  Naval  Postgraduate  School,  R.  J.  Mckee 
[Ref.  23]  confirmed  the  presence  of  Taylor-Gortler  vortices 
in  a  curved  rectangular  channel  of  aspect  ratio  40.   M. 
Durao  [Ref.  24]  modified  R.  J.  Mckee's  channel  and  modeled 
it  as  infinite  parallel  plates  with  the  outer  wall  heated 
and  inner  wall  adiabatic.   M.  Durao  and  J.  Ballard  [Refs. 
24,  25]  investigated  the  effects  of  Taylor-Gortler  vortices 
on  the  heat  transfer  in  a  straight  and  curved  test  section 
for  laminar  flows.   R.  Holihan,  Jr.  [Ref.  26]  reported 
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results  for  laminar  and  transition  flows.   S.  Daughety 
[Ref.  27]  reported  results  for  turbulent  flows  and  J.  Wilson 
[Ref.  28]  reported  results  for  transition  and  turbulent 
flows.   Each  time  Mckee's  apparatus  was  used  it  was  modified 
slightly  using  prior  recommendations  to  improve  on  the 
accuracy  of  the  results. 
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II.   INTENT  OF  STUDY 

The  purpose  of  this  investigation  was  to  build  an 
improved  curved  rectangular  duct  of  large  aspect  ratio  which 
could  be  configured  with  a  variety  of  boundary  conditions, 
and  yield  more  accurate  data  then  previously  attained.   By 
combining  new  ideas  and  experimental  equipment  with  past 
experiences  and  recommendations  of  R.  Mckee  and  others 
previously  mentioned  [Refs.  23,  24,  25,  26,  27,  28],  the 
effect  of  streamwise  curvature  on  the  heat  transfer  rate  in 
a  curved  rectangular  channel  could  be  studied  in  greater 
detail  while  providing  a  smaller  uncertainty  in  the  results. 

Some  of  the  improvements  incorporated  in  the  new  channel 
include:   the  ability  to  heat  one  or  both  walls  over  a  wide 
range  of  constant  heat  fluxes;  the  ability  to  heat  both 
walls  at  different  constant  heat  fluxes;  the  use  of  more 
calibrated  thermocouples  which  will  allow  greater  precision 
in  reading  surface  and  air  stream  temperatures,  combined 
with  a  better  arrangement  which  will  allow  the  calculation 
of  the  local  heat  transfer  coefficient  and  Nusselt  number; 
a  longer  transition  region  at  the  outlet  which  will  allow 
for  more  mixing  and  a  more  accurate  outlet  temperature 
measurement.   The  use  of  an  outer  and  inner  section  allowed 
both  wall  surfaces  to  be  specially  constructed  with  smooth 
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transitions  between  wall  and  heater  surfaces  in  both 
straight  and  curved  test  sections,  preventing  unexpected 
tripping  from  laminar  to  turbulent  flow.   Finally,  by 
constructing  the  channel  with  materials  having  low  thermal 
conductivities  and  avoiding  any  heat  sinks  along  the  channel 
walls,  heat  losses  could  be  minimized. 

Taylor-Gortler  vortices  are  known  to  enhance  heat 
transfer  in  curved  channels  for  transition  and  turbulent 
flows.   This  enhancement  in  the  transfer  of  heat  has  been 
attributed  to  the  secondary  flow  velocity  components  of  the 
Taylor-Gortler  vortices.   These  secondary  components 
transport  the  heated  fluid  from  the  outer  concave  wall, 
inward  toward  the  opposite  wall  of  the  channel,  displacing 
the  cooler  fluid  particles  and  causing  them  to  move  toward 
the  heated  concave  wall.   It  was  expected  that  similar 
results  would  be  observed  in  this  study. 

This  investigation  was  conducted  using  a  single  channel 
with  a  rectangular  cross-section  and  constant  aspect  ratio, 
again  defining  aspect  ratio  as  the  spanwise  distance  divided 
by  the  channel  height.   The  channel  incorporates  a  straight 
and  a  curved  test  section.   Each  test  section  consists  of 
opposing  walls  which  can  be  independently  heated.   The 
results  obtained  in  each  test  section  at  approximately  the 
same  flow  rates  were  compared  in  an  effort  to  determine  the 
effects  of  the  Taylor-Gortler  vortices  on  the  heat  transfer 
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rate.   Also,  the  results  of  both  the  straight  and  curved 
outer  sections  were  compared  to  the  results  of  Daughety 
[Ref.  27],  and  Wilson  [Ref.  28]  to  compare  the  new  channel 
to  the  old  one. 

The  straight  section  results  of  this  study  were  compared 
to  the  Di ttus-Boelter  correlation  using  hydraulic  diameter 
[Ref.  29],  and  the  results  of  Kays  and  Leung  for  turbulent 
and  transition  flow  in  annular  passages  [Ref.  21]. 

The  curved  section  results  of  this  study  were  compared 
with  the  results  of  Brinich  and  Graham  [Ref.  22]  for 
turbulent  flow  in  a  rectangular  curved  channel,  and  with 
Kreith  [Ref.  16]  for  turbulent  flows  in  concave  and  convex 
curved  channels. 
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III.   EXPERIMENTAL  SETUP 

A.   DESCRIPTION  OF  THE  APPARATUS 

To  meet  the  objectives  of  this  investigation,  a  channel 
of  rectangular  cross-section  and  constant  aspect  ratio  was 
constructed.   The  rectangular  channel  has  four  electrically 
heated  wall  sections,  composed  of  a  rubber  heater  bonded  to 
a  thin  copper  plate.   All  four  heated  wall  sections  are 
30.48  centimeters  long  and  25.4  centimeters  wide,  with  a 
heated  area  of  774.2  square  centimeters.   The  rectangular 
channel  is  0.635  centimeters  high  and  25.4  centimeters  wide 
resulting  in  an  aspect  ratio  of  40.   The  cross-sectional 
area  of  the  channel  is  16.13  square  centimeters.   The  wetted 
perimeter  is  52.07  centimeters,  and  the  hydraulic  diameter 
is  1.239  centimeters.   Figure  3  shows  a  cross-sectional  view 
presenting  both  the  straight  and  curved  sections  of  the 
channel . 

The  channel,  shown  in  Figure  4,  has  a  straight  entrance 
region  of  76.20  centimeters,  which  insures  hydrodynamically 
fully  developed  flow  prior  to  entering  the  straight  heated 
test  sections.   The  entrance  region  is  followed  by  two 
parallel  30.48  centimeter  straight  heated  test  sections, 
which  oppose  each  other,  one  in  each  wall.   Following  the 
straight  heated  test  sections  is  a  15.32  centimeter  straight 


25 


p 


(A 


;.;.;. ;. ;.-.-: 


CO 

CD 

C 

CO 

I 

o 

CO 

■i^B 

k. 

+* 

O) 

c 

CD 

CO 

■«■ 

co 

X) 

— 

X 

X 

_Q 

3 

0) 

0) 

3 

03 

DC         i 


4-rf 

0 

CO 

O 

CO 

Q. 

a 

k. 

CO 

0 

a 

CO 

a 

+■» 

o 

CD 

O 

2 

CO 

c 
c 
id 

u 

4-1 

o 
a) 

•H 

> 


C 

o 

•H 
+J 

o 

<D 
0) 
I 

(fl 
W 

o 

5-1 
U 

0 
■H 
+J 

fC 

e 

CO 
-C 

o 

in 


co 
CO 

CP 
•H 

Pn 


26 


C 

c 
<c 
,C 
u 

u 

03 

r-l 

cn 

c 

(T3 
-U 
O 
<D 

05 

U-l 

O 
U 

•H 
4-) 

I 
<u 

u 


cn 

•H 


27 


exit  region,  that  completes  the  straight  test  section.   The 
straight  test  section  is  immediately  followed  by  the  curved 
test  section  which  subtends  a  180  degree  arc,  with  the  inner 
wall  having  a  radius  of  curvature  of  29.7  centimeters.   The 
curved  test  section  has  a  52.63  centimeter  unheated  convex 
curved  section  prior  to  the  curved  heated  sections.   The  two 
concentric  curved  heated  sections,  inner  one  convex  and 
outer  one  concave,  are  30.48  centimeters  long  and  the  outer 
heated  wall  subtends  an  arc  of  58.8  degrees.   The  curved 
heated  sections  are  followed  by  a  10.16  centimeter  convex 
curved  exit  region  which  completes  the  curved  test  section. 
Finally,  the  channel  is  completed  by  a  92.0  centimeter 
straight  exit  region  which  allows  sufficient  mixing  prior  to 
reading  the  outlet  temperature.   The  entire  channel  was  then 
thermally  insulated  with  a  plastic  foam  material. 

An  entrance  bell  constructed  of  Plexiglas  is  connected 
to  the  inlet  of  the  straight  entrance  section  of  the 
channel.   It  was  designed  and  manufactured  according  to  ASME 
nozzle  standards,  with  an  elliptical  curved  base  on  a  major 
axis  equal  to  25.4  centimeters  and  a  minor  axis  of  2.54 
centimeters.   A  fiberglass  filter  was  attached  to  the 
entrance  nozzle  to  prevent  foreign  matter  from  entering  the 
channel.   An  aluminum  exhaust  nozzle  was  attached  to  the 
exit  of  the  channel  and  directed  the  flow  from  the  channel 
into  a  two  inch  diameter  pvc  pipe.   The  piping  contained  a 
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mass  flow  rate  measuring  section  which  could  be  fitted  with 
standard  ASME  concentric  orifices. 

A  total  of  one  hundred  and  eighty-nine  copper-constantan , 
30  gauge  thermocouples  are  located  throughout  the  test 
apparatus.   Only  ninety-nine  thermocouples  are  read  at 
anyone  time,  through  seventy-five  data  channels,  to  record 
the  desired  temperatures.   Each  heated  wall  section  contains 
thirty  thermocouples  which  are  read  individually  and  then 
averaged  to  give  the  wall  section  surface  temperature. 
Four  thermocouples  span  the  width  of  the  channel  at  the 
entrance  and  exit  to  give  the  average  entrance  and  exit  air 
temperatures.   Three  groups  of  five  thermocouples  are 
connected  in  parallel  and  inserted  between  each  layer  of 
insulation  at  each  of  the  heated  test  sections.   Each  group 
of  five  thermocouples  only  uses  one  data  channel,  and  these 
extra  thermocouples  comprise  the  difference  between  the 
number  of  thermocouples  read  and  the  data  channels 
available.   Finally,  there  is  one  thermcouple  to  measure  the 
air  temperature  at  the  orifice. 

B.   DETAILS  OF  THE  CONSTRUCTION 

Each  of  the  four  heated  wall  sections  was  composed  of  a 
copper  plate  laminated  to  a  wire  wound  silicone  rubber 
heater.   The  rubber  heater  was  made  by  Watlow  Industries. 
There  are  thirty  copper-constantan  thermocouples  placed 
between  the  copper  plate  and  rubber  heater.   Figure  5  shows 
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Figure  5.   Thermocouple  Placement  in  Copper  Plates 
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the  placement  of  the  thermocouples  on  the  back  of  the  copper 
plate.   The  copper  plate  is  30.48  centimeters  long,  26.72 
centimeters  wide,  and  0.155  centimeters  thick,  but  only  the 
middle  25.4  centimeters  of  copper  plate  width  is  heated  and 
centered  in  the  channel.   The  back  of  the  copper  plate  has  a 
pattern  of  slots  milled  0.076  centimeters  deep.   The  pattern 
of  slots  for  the  two  inner  plates  is  different  from  the 
pattern  of  slots  for  the  two  outer  plates.   See  Figures  6, 
7,  8,  and  9.   The  calibrated  thermocouples  were  numbered  and 
expoxied  into  the  slots,  as  shown  in  Figures  10  and  11.   The 
rubber  heater  has  the  same  length  and  width  dimensions  as 
the  copper  plate,  but  only  heats  the  middle  25.4  centimeters 
of  width.   This  results  in  a  heated  wall  area  of  774.2 
square  centimeters.   The  excess  copper  plate  is  sandwiched 
between  the  wall  and  the  side  spacer  thereby  helping  to 
support  the  copper  plate  in  the  channel.   Dow  Corning  3140 
RTV  was  used  to  adhere  the  rubber  heater  to  the  back  of  the 
copper  plate,  and  complete  the  heated  section.   Figure  12 
and  13  represent  a  completed  test  wall. 

The  channel  was  built  on  a  carefully  constructed  wood 
frame  with  steel  rails,  that  also  was  used  to  support  the 
entire  channel.   The  inner  wall  was  composed  of  two  0.343 
centimeter  thick  and  29.85  centimeter  wide  pieces  of 
Plexiglas  laminated  together,  that  had  spaces  milled  out  for 
the  heated  copper  sections.   The  inner  wall  was  constructed 
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Figure  6.   Photograph  of  Milled  Thermocouple  Slots 


Figure  7.   Photograph  of  Thermocouple  Placement  Holes 
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Figure  8.   Photograph  of  Both  Inner  Plates  Showing 
Thermocouple  Slot  Pattern. 
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Thermocouple  Slot  Pattern. 
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Figure  12.   Photograph  of  Convex  Surface  of  Copper  Plate 


Figure  13.   Photograph  of  Completed  Heated  Test  Section 
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by  heating  and  forming  the  Plexiglas  around  the  wood  frame. 
The  Plexiglas  wall  was  secured  to  the  frame  by  screws  tapped 
into  the  steel  rails.   The  heated  copper  sections  were  then 
fitted  into  the  machined  spaces  and  two  0.635  centimeter 
thick  and  2.223  centimeter  wide  Plexiglas  spacers,  that 
served  as  the  sides  of  the  channel  and  support  for  the 
copper  plates,  were  attached  by  screws  to  the  inner  wall  and 
steel  rails.   The  two  spacers  were  25.4  centimeters  apart 
and  formed  the  width  of  the  channel.   Additionally,  five 
flat  head  screws  spanning  the  width  of  the  channel  at  the 
joints  where  the  copper  and  Plexiglas  met  were  used  to 
secure  the  copper  plates  and  maintain  an  even  joint.   The 
seams  were  then  worked  to  a  smooth  transition  by  applying 
"Bondo"  to  the  joints  and  sanding.   See  Figures  14,  15,  16, 
and  17. 

The  outer  wall  was  constructed  in  a  similar  manner,  but 
using  0.318  centimeter  thick  and  29.85  centimeter  wide 
pieces  of  Lexan.   Lexan  was  used  on  the  outer  wall  because 
of  its  flexibility  and  ease  of  bending  without  heat,  which 
would  be  required  when  attaching  the  outer  wall  to  the  inner 
wall  and  frame.   Also,  the  outer  wall  was  formed  in  an  outer 
wood  frame  so  that  the  transition  joints  of  the  outer  wall 
could  be  worked  smooth  like  the  inner  wall.   See  Figures  18 
and  19.   Both  walls  were  constructed  so  that  no  seams  or 
butt  joints  existed  on  the  inside  of  the  channel,  except 
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Figure  14.   Photograph  of  Channel  Frame 


Figure  15. 


Photograph  of  Bottom  View  of  Inner  Straight 
Test  Section. 
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Figure  16.   Photograph  of  Inner  Straight  Wall  at 
Copper-Plexiglas  Joint. 


Figure  17. 


Photograph  of  Inner  Curved  Wall  at 
Copper-Plexiglas  Joint. 
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Figure  18. 


Photograph  of  Outer  Frame  and  Lexan  Wall 
Showing  Straight  Test  Section. 
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Figure  19.   Photograph  of  Outer  Curved  Test  Section 
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where  the  heated  copper  section  joined  the  Plexiglas  and 
Lexan  walls,  and  all  eight  of  those  joints  had  a  smooth 
transition. 

Once  both  walls  were  completed,  the  outer  wall  and 
support  frame  was  positioned  onto  the  inner  wall  and  main 
frame.   The  outer  wall  was  aligned  onto  the  inner  wall  and 
the  outer  wall  support  frame  removed.   The  outer  wall  was 
then  bolted  to  the  inner  wall  and  frame  with  steel  bars 
running  the  length  of  the  channel  used  as  a  spacer  to  insure 
even  compression.   Once  this  was  completed  the  inlet  and 
outlet  nozzles  were  attached  to  the  channel. 

The  sides  of  the  channel  and  all  connections  between  the 
channel,  the  pvc  piping,  and  the  orifice  were  sealed  with 
Macklanburg-Duncan  Silicone  Rubber  Sealant  Caulk  to  ensure 
there  was  no  leakage  of  air  into  the  channel  or  the  piping 
which  could  effect  the  recorded  temperatures  or  mass  flow 
rate.   The  entire  channel  was  thermally  insulated  with  three 
layers  of  1/2  inch  Armstrong  Armaflex  22  Sheet  Insulation, 
a  flexible  foamed  plastic  material.   By  positioning 
thermocouples  between  these  insulation  layers,  the  heat  lost 
through  the  insulation  could  be  computed  at  each  test 
section.   The  heat  loss  through  insulation  surrounding  the 
rest  of  the  channel  due  to  the  heated  air  flowing  through 
the  channel  was  a  small  fraction  of  the  total  power  supplied 
to  the  heater.   The  insulation  was  held  in  place  by  the  use 
of  silver  duct  tape  and  velcro  strips  attached  to  the  frame. 
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C.   INSTRUMENTATION 

The  data  acquisition  system  used  for  this  experiment  was 
a  Hewlett  Packard  3054A  Automatic  Data  Acquisition/Control 
System  consisting  of  a  3456A  Digital  Voltmeter,  a  3497A  Data 
Acquisition/Control  Unit,  and  four  input  panels  capable  of 
reading  any  mix  of  eight  voltages  or  copper-constantan 
thermocouples.   Also  used  in  conjunction  with  this  data 
acquisition  system,  were  a  Hewlett  Packard  9826  Computer  and 
2671A  Printer.   The  program  used  to  reduce  the  data  is 
listed  in  Appendix  A. 

The  copper-constantan  thermocouples  were  calibrated 
against  a  platinum  resistance  thermometer  using  a  ROSEMONT 
Commutating  Bridge  model  920A,  and  a  ROSEMONT  Constant 
Temperature  Bath  model  913A.   All  thermocouples  were 
calibrated  between  20  and  85  degrees  Celsius,  and  a  second 
order  polynomial  was  curve  fitted  to  correct  each 
thermocouple. 

The  mass  flow  rate  in  the  channel  was  measured  with  a 
thin  plate  orifice  measuring  device  which  was  placed  in  the 
2  inch  pvc  piping  between  the  channel  and  the  compressor. 
In  this  study  an  orifice  with  a  diameter  of  2.7318 
centimeters  was  used  to  measure  hydraulic  Reynold's  numbers 
above  5000.   The  orifice  device  was  installed  in  accordance 
with  ASME  Power  Test  Code  [Ref.  30],  with  3/8  inch  pressure 
taps  placed  on  either  side  of  the  orifice  at  1  and  1/2 
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diameters,  and  one  thermocouple  inserted  in  the  two  inch  pvc 
piping  down  stream  of  the  orifice  to  record  the  temperature 
of  the  air  flowing  through  the  orifice.   Pressure  measure- 
ments were  taken  with  Meriam  vertical  manometers,  one  whose 
fluid  was  water  with  a  0  to  60  inch  range,  another  using 
water  with  a  0  to  30  inch  range,  and  a  third  whose  fluid 
was  mercury,  calibrated  to  read  inches  of  water  with  a  range 
of  0  to  415  inches.   Atmospheric  pressure  was  measured  with 
a  Princo  precision  barometer.   The  pressure  readings  were 
used  to  calculate  the  mass  flow  rate  of  the  working  fluid  as 
outlined  in  [Refs.  30,  31]. 

Power  was  supplied  to  the  heated  test  section  by  using 
two  Lambda  Model  LK345A  FM  60  volt  power  supplies  in  series. 
The  power  supplies  are  rated  at  6.0  amps,  maximum  at  40 
degrees  Celsius.   Since  the  electrical  resistance  of  the 
rubber  heater  is  not  constant,  but  varies  slowly  with 
temperature,  a  precision  resistor  was  calibrated  and 
connected  in  series  with  the  heater  to  allow  the  calculation 
of  the  instantaneous  power  being  supplied.   Both  the  outer 
and  inner  heated  test  sections  have  their  own  circuitry 
and  precision  resistor.   The  value  of  the  precision 
resistor  for  the  outer  test  sections  is  2.0173  ohms,  and 
the  value  for  the  inner  test  sections  is  2.0084  ohms.   See 
Figure  20  for  diagram  of  circuit  used  to  measure 
instantaneous  power. 
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Figure  20.   Power  Measurement  Circuit 
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The  working  fluid  was  air  at  room  temperature,  which  was 
drawn  through  the  channel  by  an  electrically  driven  Spencer 
Turbo  Compressor,  rated  at  30  horsepower  at  3500  rpm,  and 
550  cubic  feet  per  minute  at  70  degrees  fahrenheit  and  one 
atmosphere.   Photographs  of  the  channel  and  associated  test 
equipment  are  shown  in  Figures  21,  22,  23,  and  24. 
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IV.   EXPERIMENTAL  PROCEDURES 

The  experimental  procedures  followed  were  the  same  for 
the  straight  and  curved  test  sections.   The  procedure  for 
starting  up  the  experimental  apparatus  is  as  follows: 

(1)  turn  on  the  computer  and  data  acquisition  system 

(2)  close  globe  valve  at  channel  exit 

(3)  close  valve  to  manometer  reading  PI 

(4)  set  compressor  blast  gate  to  1/8  open 

(5)  start  turbo  compressor 

(6)  open  intake  valve  on  compressor  1/4  to  1  turn 

(7)  open  globe  valve  at  channel 

(8)  open  manometer  valve  slowly  to  ensure  PI  does  not 
exceed  manometer  scale 

(9)  turn  on  and  adjust  power  supply  to  the  heater 

(10)  use  intake  valve  on  compressor  to  control  mass  flow 
rate  through  the  channel 

It  was  necessary  to  determine  the  approximate  voltage 

setting  required  to  bring  the  heated  wall  to  approximately 

50  degrees  Celsius.   This  heated  wall  temperature  would 

ensure  a  large  temperature  difference  between  the  inlet  and 

outlet  temperatures,  and  the  heated  and  unheated  wall 

temperatures,  while  not  damaging  the  apparatus.   The  larger 

the  temperature  difference  the  smaller  the  uncertainty  in 

the  results.   The  required  power  setting  was  found  by 


51 


monitoring  the  heated  plate  temperature  locally  from  the 
data  acquisition  system  while  the  voltage  to  the  heater  was 
adjusted . 

To  compute  the  instantaneous  power  supplied  to  the 
heater  the  following  relationship  was  used: 


V   V 
H  VPR 


R 


PR 


The  precision  resistor  voltage  and  the  heater  voltage  were 
read  by  the  data  acquisition  system  and  the  resistance  of 
the  precision  resistor  was  a  known  constant.   A  diagram  of 
the  circuit  used  was  shown  in  Figure  16. 

The  mass  flow  rate  was  monitored  by  the  two  manometers 
at  the  orifice.   For  each  run,  the  atmospheric  pressure,  the 
pressure  difference  across  the  orifice,  and  the  pressure 
upstream  of  the  orifice  were  entered  manually  into  the 
computer  for  the  mass  flow  rate  calculation.   All  other  data 
was  acquired  by  the  data  acquisition  system  directly  from 
the  experimental  apparatus. 

Once  the  apparatus  was  started  properly,  sufficient  time 
was  required  to  reach  steady  state  prior  to  gathering  data. 
Preliminary  runs  were  performed  to  determine  the  time 
required  for  the  test  apparatus  to  reach  steady  state.   For 
this  experiment  at  least  two  hours  were  allowed  for  a  cold 
test  apparatus  to  reach  steady  state,  and  at  least  one  hour 
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if  the  apparatus  was  already  running.   After  waiting  the 
required  time,  data  was  taken  at  ten  minute  intervals  to 
ensure  that  a  steady  state  condition  had  been  achieved.   The 
criteria  for  steady  state  was  based  on  three  variables: 

(1)  the  difference  between  inlet  and  outlet  temperature 

(2)  the  heated  test  section  wall  temperature 

(3)  the  unaccounted  heat  loss  in  the  system 

When  the  first  two  of  these  quantities  varied  by  less  than 
two  percent  over  a  ten  minute  interval  and  the  third  was 
less  than  ten  percent  of  the  input  power,  it  was  considered 
that  a  steady  state  condition  had  been  achieved. 

Experiments  were  conducted  at  specific  mass  flow  rates 
of  air  corresponding  to  Reynolds  numbers  between  7000  and 
22000.   The  Reynolds  number  was  based  on  the  hydraulic 
diameter  of  the  channel.   For  each  mass  flow  rate,  tempera- 
ture and  voltage  data  were  taken  automatically  by  the  data 
acquisition  system  and  immediately  reduced  and  printed.   The 
variables  that  were  measured  and  used  in  this  investigation 
were  : 

(1)  the  temperature  of  the  air  entering  the  channel  (T.  ) 

(2)  the  temperature  of  the  air  leaving  the  channel  (TQUt) 

(3)  the  temperature  of  the  heated  wall  in  the  test 
section  (TWQ  or  Twl ) 

(4)  the  temperature  of  the  unheated  wall  in  the  test 
section  (Twi  or  Twq ) 

(5)  the  temperature  between  each  layer  of  insulation  at 
the  test  section  (T.   ) 

1  lib 


53 


(6)  the  temperature  of  the  air  at  the  orifice  (Torf) 

(7)  the  voltage  across  the  precision  resistor  (VpR) 

(8)  the  voltage  across  the  rubber  heater  (VH) 

(9)  the  pressure  upstream  of  the  orifice  (PI) 

(10)  the  difference  in  the  upstream  and  downstream 
pressure  across  the  orifice  (AP) 

(11)  the  atmospheric  pressure  (p^<pM^ 

The  results  obtained  from  the  straight  test  sections 
served  as  a  baseline  for  the  comparison  of  the  curved 
section  results.   An  effort  was  made  to  match  the  mass  flow 
rates  for  each  test  section  as  closely  as  possible. 
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V.   PRESENTATION  OF  DATA 

A.   DATA  REDUCTION 

The  constant  heat  flux  surface  for  the  heated  wall  of 
each  test  section  was  provided  by  the  uniform  electrical 
resistivity  of  the  rubber  heater.   The  insulated  unheated 
wall  was  considered  adiabatic  since  the  heat  losses  through 
that  wall  were  negligible.   The  straight  portion  of  the 
channel  upstream  of  the  straight  heated  test  section  was  of 
sufficient  length  to  ensure  that  the  air  flow  was  hydro- 
dynamically  fully  developed  for  the  flow  velocities  of  this 
experiment.   The  straight  portion  of  the  channel  downstream 
of  the  curved  section  proved  to  be  of  sufficient  length  to 
ensure  that  the  air  flow  exiting  the  channel  was  thoroughly 
mixed  and  that  the  measured  outlet  air  temperature  (Tout) 
was  an  average  bulk  temperature,  as  long  as  the  hydraulic 
Reynolds  number  of  the  channel  was  less  than  17000. 

The  analysis  of  the  data  required  several  quantities  to 
be  defined.   The  heat  convected  to  the  air  was  calculated 
using  the  equation: 

Qair  =  mCpair  (Tout  "  Tin) 

where  'C   ,  '  is  the  specific  heat  of  the  air  at  constant 
pa  1  r  r 

pressure  and  'm'  is  the  mass  flow  rate  of  the  air. 
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The  average  heat  transfer  coefficient  between  the  heated 
wall  and  the  flow  of  air  in  the  channel  was  defined  by  the 
equation: 


h   = 


Q  • 
air 


where  'Q  ,  '  is  defined  above,  'AD, '  is  the  area  of  the 

air  rL 

heated  copper  plate  in  the  test  section,  and  'AT'  is  the 
difference  between  the  arithmetic  average  heated  wall 
temperature  (T  )  and  the  average  bulk  air  temperature 
(T,  , k) .   The  average  bulk  air  temperature  is  defined  as  the 
arithmetic  mean  of  the  entrance  and  exit  air  temperatures 

(Tin'  Tout>- 

The  average  Nusselt  number  was  then  calculated  using: 


h  D 


Nu 


hd 


hd 


K 


air 


In  this  equation  'D,  ,'  is  the  hydraulic  diameter  and  'K  .  ' 
^         nd  air 

is  the  thermal  conductivity  of  the  air. 

The  hydraulic  Reynolds  number  was  calculated  for  each 
test  run  as  follows: 


m  D 


Re 


hd 


hd 


A   u  . 
c   air 


where  again  'm'  and  'D.  , '  are  the  mass  flow  rate  and 
hydraulic  diameter  of  the  channel,  'A  *  is  the  cross- 
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sectional  area  of  the  channel,  and  'u  .  '  is  the  dynamic 

air 

viscosity  of  the  air. 

For  the  curved  sections  the  Dean  number  was  defined  as 


m  Dc/2 

a    r~~ 

c   air 


Re 


where  'Re,'  is  the  Reynolds  number  based  on  the  channel 
height,  ' Dc '  *s  tne  channel  height,  and  'Ri '  is  the  radius 
of  curvature  of  the  inner  convex  wall. 

A  sketch  of  the  control  volume  and  a  set  of  sample 
calculations  for  one  test  run  of  the  curved  section  are 
given  in  Appendix  B.   A  sample  calculation  for  the 
uncertainty  analysis  is  given  in  Appendix  C. 

B.   RESULTS 

The  data  obtained  from  each  experimental  run  was  reduced 
by  the  computer  program  listed  in  Appendix  A,  utilizing  the 
expressions  described  in  Part  A  above.   The  major  parameters 
resulting  from  this  process  are  shown  in  Tables  I  through 
IV.   Tables  I  and  II  contain  the  straight  sections  results 
and  Tables  III  and  IV  contain  the  results  of  the  curved 
sections.   A  plot  of  the  average  Nusselt  number  versus 
hydraulic  Reynolds  number  is  given  in  Figures  25,  26,  27, 
28,  and  29  for  the  comparison  of  the  results.   "Error"  bands 
have  been  indicated  as  a  result  of  the  uncertainty  analysis. 


57 


TABLE  I 
SUMMARY  OF  STRAIGHT  OUTER  WALL  TEST  RESULTS 


Re,  Q  .  h  AT         Nu,  . 

hd         (W)r       (W/mC)        (C)  hd 


7709  108.50  41.99  28.52  20.07 

11953  135.26  61.57  25.49  29.51 

15646  157.92  74.43  24.45  35.65 

21455  217.53  87.42  27.36  41.95 


Reynolds  number  and  Nusselt  number  are  based  on 
hydraulic  diameter.   Results  for  Reynolds  number  of 
21455  based  on  orifice  air  temperature  instead  of 
outlet  air  temperature. 
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TABLE  II 
SUMMARY  OF  STRAIGHT  INNER  WALL  TEST  RESULTS 


Rehd 

Q  • 
(W}r 

h 
(W/mC) 

AT 
(C) 

Nuhd 

7628 

105.57 

41.52 

28.90 

19.85 

12169 

138. 36 

59.50 

27.42 

28.  50 

15788 

159.84 

70.24 

26.37 

33.61 

21312 

204.42 

81.82 

27.84 

39.17 

Reynolds  number  and  Nusselt  number  are  based  on 
hydraulic  diameter.   Results  for  Reynolds  number  of 
21312  based  on  orifice  air  temperature  instead  of 
outlet  air  temperature. 
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TABLE  III 
SUMMARY  OF  CURVED  CONCAVE  WALL  TEST  RESULTS 


Rehd 

De 

Q  • 

(6)r 

h 
(W/mC) 

T 
(C) 

Nuhd 

7782 

206 

115.88 

52.43 

25.51 

25.00 

11726 

311 

157.45 

68.16 

29.91 

32.53 

15813 

419 

158.17 

83.  13 

22.34 

39.78 

21246 

563 

231.31 

93.02 

27.20 

44.52 

Reynolds  number  and  Nusselt  number  are  based  on 
hydraulic  diameter.   Results  for  Reynolds  number  of 
21246  based  on  orifice  air  temperature  instead  of 
outlet  air  temperature. 
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TABLE  IV 
SUMMARY  OF  CURVED  CONVEX  WALL  TEST  RESULTS 


Rehd 

De 

Q  • 

(wfr 

h 
(W/mC) 

AT 
(C) 

Nuhd 

7757 

206 

82.37 

38.33 

24.47 

18.37 

12021 

318 

115. 50 

57.62 

23.  92 

27.65 

15908 

421 

142.54 

68.57 

24.34 

32.94 

21603 

572 

164.08 

79.28 

23.43 

38.  11 

Reynolds  number  and  Nusselt  number  are  based  on 
hydraulic  diameter.   Results  for  Reynolds  number  of 
21603  based  on  orifice  air  temperature  instead  of 
outlet  air  temperature. 
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As  would  be  expected  the  results  indicate  an  increase  in 
the  rate  of  heat  transfer  with  increasing  Reynolds  number 
for  both  the  straight  and  curved  test  sections.   In 
addition,  the  heat  transfer  rate  was  highest  in  the  concave 
curved  section  for  each  Reynolds  number  investigated,  and 
lowest  in  the  convex  curved  section. 

In  the  low  end  of  the  turbulent  regime,  at  a  hydraulic 
Reynolds  number  of  7780,  the  rate  of  heat  transfer  was  about 
twenty-two  percent  higher  in  the  concave  curved  section  than 
the  outer  straight  section.   As  the  hydraulic  Reynolds 
number  continued  to  increase,  between  11700  and  15800,  the 
rate  of  heat  transfer  was  approximately  eleven  percent 
higher  between  the  concave  curved  section  and  the  outer 
straight  section.   At  a  hydraulic  Reynolds  number  of  21000 
the  difference  in  heat  transfer  rate  decreased  to  about 
seven  percent. 

The  results  obtained  at  a  hydraulic  Reynolds  number  of 
21000  have  a  much  larger  uncertainty,  because  the  orifice 
air  temperature  was  higher  than  the  average  outlet  air 
temperature.   This  signifies  that  the  air  was  not  fully 
mixed  at  the  outlet,  because  no  heat  is  added  to  the  air 
after  the  outlet  air  temperature  is  read,  yet  the  air 
temperature  increased.   For  hydraulic  Reynolds  numbers 
greater  than  about  17000  the  channel  outlet  section  is  not 
long  enough  to  allow  the  air  to  mix  properly. 
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The  values  obtained  for  the  hydraulic  Reynolds  number  of 
approximately  21000  are  based  on  the  orifice  air  temperature 
vice  the  outlet  air  temperature.   This  resulted  in  a  Nusselt 
number  higher  than  obtained  by  using  the  the  outlet  air 
temperature,  but  still  lower  than  the  actual  Nusselt  number, 
because  of  the  small  heat  loss  between  the  outlet  nozzle  and 
the  thermocouple  at  the  orifice.   This  heat  loss  cannot  be 
accounted  for  and  the  temperature  used  as  an  outlet  air 
temperature  was  still  lower  than  the  actual  outlet  air 
temperature . 

In  comparing  the  concave  curved  wall  to  the  convex 
curved  wall  the  rate  of  heat  transfer  was  about  thirty-six 
percent  higher  in  the  concave  section  for  a  hydraulic 
Reynolds  number  of  7780.   For  hydraulic  Reynolds  numbers 
between  11700  and  15800  the  rate  of  heat  transfer  was  twenty 
percent  higher  in  the  concave  section,  and  at  a  hydraulic 
Reynolds  number  of  21200  the  rate  of  heat  transfer  was 
seventeen  percent  higher. 
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VI .   DISCUSSION  AND  CONCLUSIONS 

J.  Wilson  [Ref.  28],  in  his  study,  reported  that  the 
actual  presence  of  Taylor-Gortler  vortices  were  verified  by 
the  observation  of  the  liquid  crystal  isotherm  distribution 
in  the  curved  test  section.   He  further  reported,  that  at 
low  values  of  hydraulic  Reynolds  number  and  through  most  of 
the  transition  regime,  liquid  crystal  bands  were  oriented  in 
the  streamwise  direction.   As  the  Reynolds  number  increased 
turbulent  forces  began  to  dominate  in  the  curved  section  and 
less  evidence  could  be  detected  for  the  existence  of  the 
vortices . 

The  high  aspect  ratio  of  the  channel  provided  the  basis 
for  the  assumption  that  the  experimental  apparatus  should 
behave  like  infinite  parallel  plates.   Holihan's  data  [Ref. 
26]  tended  to  verify  this  assumption.   His  experimental  data 
in  the  laminar  flow  region,  approached  the  theoretical  limit 
for  average  Nusselt  number  for  parallel  plates  with  one  wall 
heated  at  a  constant  heat  flux  and  the  opposite  wall 
adiabatic  [Ref.  32],   Based  on  the  assumptions  mentioned 
above,  comparisons  were  made  with  other  analytical  and 
experimental  results  that  were  of  the  same  or  similar 
problems. 

A  comparison  of  the  experimental  results  of  this  study 
and  the  experimental  results  of  Daughety  [Ref.  27]  and 
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Wilson  [Ref.  28]  are  shown  in  Figure  30  for  the  curved 
section,  and  Figure  31  for  the  straight  section.   The 
channel  aspect  ratio  was  40  for  each  of  these  studies  and 
the  experimental  procedures  were  similar,  with  the  exception 
of  a  slight  modification  in  Wilson's  procedures.   Wilson 
[Ref.  28]  made  a  correction  to  determine  the  temperature 
drop  across  the  liquid  crystal  sheet  covering  the  heated 
wall.   By  applying  this  temperature  drop  to  the  measured 
temperature,  the  true  temperature  of  the  heated  wall  would 
be  known.   After  reviewing  Wilson's  data  and  experimental 
procedure  it  was  detected  that  Wilson  was  over  estimating 
the  temperature  drop  across  the  liquid  crystal  by 
approximately  five  to  seven  percent,  and  this  lead  to  an 
over  estimation  of   T  of  ten  to  twelve  percent.   This 
resulted  in  an  over  estimation  of  the  average  heat  transfer 
coefficient  and  the  average  Nusselt  number  of  between  ten 
and  twelve  percent.   When  Wilson's  data  is  lowered  by  ten 
percent,  both  his  data  and  Daughety's  data  fall  within  the 
uncertainty  bands  of  the  present  study  as  shown  in  Figures 
32  and  33.   Differences  were  attributed  to  the  inability  to 
acccount  for  all  of  the  heat  transfer  processes  and  the 
limitations  inherent  in  any  experimental  work. 

Ballard  [Ref.  25]  reported  an  increase  in  the  rate  of 
heat  transfer  between  a  concave  wall  and  a  straight  outer 
wall  of  eleven  percent  for  laminar  flows.   Holihan  [Ref.  26] 
reported  an  increase  of  fifteen  percent  for  laminar  flows 
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and  thirty  percent  for  transition  flows.   Daughety  [Ref. 
27],  reported  an  increase  of  twenty  percent  in  the  heat 
transfer  rate  for  turbulent  flows.   Wilson  [Ref.  28], 
reported  an  increase  of  forty-five  percent  for  laminar 
flows,  twenty  percent  for  transition  flows,  thirteen  percent 
for  hydraulic  Reynolds  numbers  between  5400  and  12000,  and 
about  twenty-six  percent  at  a  hydraulic  Reynolds  number  of 
24000.   This  compares  very  favorably  with  the  present 
results,  except  for  the  data  at  a  hydraulic  Reynolds  number 
of  21000.   No  conclusion  can  be  made  about  the  present 
results  at  a  hydraulic  Reynolds  number  of  21000  because  of 
the  high  uncertainty  of  the  data  as  previously  mentioned. 
This  area  needs  further  experimentation  before  accurate 
conclusions  can  be  made  about  the  heat  transfer  rate  at  high 
hydraulic  Reynolds  numbers. 

Kreith  [Ref.  16]  reported  an  increase  in  the  rate  of 
heat  transfer  along  a  concave  wall  as  compared  to  a  convex 

wall  of  twenty-five  to  sixty  percent  for  hydraulic  Reynolds 

4       6 
numbers  between  10   to  10  .   This  also  compares  favorably 

with  the  present  results  which  indicate  an  increase  of 

seventeen  to  thirty-eight  percent  when  comparing  concave 

wall  to  convex  wall. 

A  comparison  of  analytical  and  experimental  results  with 

those  of  the  present  study  for  the  straight  test  section  is 

shown  graphically  in  Figure  34.   The  Di ttus-Boel ter  equation 
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[Ref.  29]  for  heat  transfer  in  a  straight  tube  for  constant 
wall  temperature  is  plotted  from: 


Nu  =  0.02Re 


0.8 


For  a  Prandtl  number  equal  to  0.71.   Additionally,  the 
results  of  Kays  and  Leung  [Ref.  21]  for  heat  transfer  in 
annular  passages  are  shown  for  r*  =  1.0,  which  equates  to 
parallel  plates.   When  examined  these  results  compare 
favorably  with  the  present  data. 

The  curved  section  data  is  plotted  in  Figure  35,  and  is 
compared  to  the  experimental  work  of  Brinich  and  Graham 
[Ref.  22].   The  accuracy  of  the  data  points  used  from 
Brinich  and  Graham  are  subject  to  some  error,  in  that  the 
actual  values  are  not  given  in  their  study  and  had  to  be 
taken  from  a  plot  of  Stanton  number  versus  Reynolds  number. 

After  carefully  examining  the  past  and  present  data  the 
trends  and  results  are  consistent.   Heat  transfer  rates 
increase  with  increasing  Reynolds  numbers.   Also,  the 
concave  curved  section  enhances  the  rate  of  heat  transfer 
over  the  convex  curved  section  and  the  straight  test 
sections  for  laminar,  transition,  and  turbulent  flows.   The 
presence  of  Taylor-Gortler  vortices  in  the  laminar  and 
transition  regime  seems  to  be  the  major  contributing  factor 
to  this  enhancement.   Whether  the  Taylor-Gortler  vortices 
continue  to  exist  in  the  turbulent  flow  regime  remains  to  be 
seen . 
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VII.   RECOMMENDATIONS 

The  amount  of  experimental  work  remaining  in  the  study 
of  heat  transfer  in  curved  channels  is  still  significant. 
The  present  experimental  apparatus  has  been  tested  in  the 
turbulent  regime,  and  the  test  data  proves  that  accurate 
results  can  be  obtained  for  hydraulic  Reynolds  numbers  below 
17000,  slightly  higher  when  heating  the  straight  test 
section.   Additional  experiments  into  the  laminar, 
transition,  and  turbulent  region  should  be  conducted  to 
better  understand  the  heat  transfer  processes  occurring  in 
the  straight  and  curved  sections.   By  remaining  below  a 
hydraulic  Reynolds  number  of  17000  accurate  data  can  be 
obtained  and  significant  conclusions  can  be  drawn  in  the 
comparison  of  the  curved  sections  to  the  straight  section, 
and  the  concave  curved  wall  to  the  convex  curved  wall. 

To  improve  the  accuracy  of  the  data  and  results  some 
slight  modifications  to  the  present  channel  should  be  done. 
First,  the  straight  outlet  section  of  the  channel  should  be 
lengthened  to  ensure  proper  mixing  of  the  heated  outlet  air. 
Second,  some  sort  of  honey  combing,  screening,  or  other  type 
of  device  should  be  placed  in  the  outlet  section  to  enhance 
the  mixing  of  the  heated  outlet  air.   Finally,  extra 
thermocouples  should  be  placed  at  the  outlet  and  read  more 
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then  once  during  a  data  run  to  decrease  the  uncertainty  of 

the  outlet  air  temperature  (T   .)  and  the  uncertainty  of  the 

r  out  J 

difference  between  the  inlet  and  outlet  air  temperature.   By 
applying  these  recommendations  the  uncertainty  of  the 
Nusselt  number  can  be  reduced  to  four  percent  or  less. 

Further  experiments  at  hydraulic  Reynolds  numbers 
greater  than  17000  also  need  to  be  examined.   This  can  be 
accomplished  if  the  above  recommendations  are  taken  along 
with  one  other.   The  channel  walls  need  outer  support  to 
prevent  the  center  of  the  channel  from  collapsing  at  the 
high  hydraulic  Reynolds  numbers.   Slight  collapsing  was 
noticed  at  hydraulic  Reynolds  numbers  greater  than  20000. 
This  collapsing  was  caused  by  the  laminating  of  two  0.3175 
centimeter  thick  walls  which  were  not  as  rigid  as  the  single 
0.635  centimeter  thick  walls  used  previously.   The  actual 
hydraulic  Reynolds  numbers  that  can  be  attained  and  still 
yield  accurate  results  will  be  limited  by  the  power  of  the 
compressor  and  the  rigidity  of  the  channel's  walls. 

Finally,  once  a  strong  data  base  can  be  obtained  for 
each  test  section  further  experimentation  can  be  conducted 
with  different  boundary  conditions.   For  example,  both 
curved  walls  heated  at  the  same  constant  heat  flux.   Once 
this  can  be  done  a  true  comparison  of  the  heat  transfer  rate 
of  the  total  curved  test  section  can  be  compared  to  the 
total  straight  test  section,  i.e.,  both  curved  walls  versus 
both  straight  walls. 
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APPENDIX    A:       COMPUTER    PROGRAM 


1000  !      FILE   NAME:    DRP2A 

1010  !      LAST   REVISION   DATE:      21    OCTOBER    198S 

1020  !      AUTHOR:      LT.    GEORGE   G.    GALYO 

1030! 

1040  !   PROGRAM  FOR  GATHERING  AND  REDUCING  DATA  WHEN  HEATING  ONE  COPPER  PLATE, 

1050  !   WITH  10  MINUTES  BETWEEN  RUNS.   PROGRAM  IS  MODIFIED  TO  ACCOUNT  FOR 

1060  !   BROKEN  THERMOCOUPLE  CHANNELS  5,  16,  21,  AND  26  WHICH  EFFECT  OUTER 

1070  !   CURVED  HEATER,  AND  THERMOCOUPLE  CHANNEL  47  WHICH  EFFECTS  INNER  CURVED 

1080  !   HEATER. 

1090! 

1100  COM  /Co/  D(7),Aa(76,2) 

1110  DIM  Emf(79),T(79) 

1120! 

1130  !   CORRELATION  FACTORS  FOR  EMF  TO  DEGREES  CELSIUS  FOR  CU-CO  T/C. 

1140! 

1150  DATA   0.10086091,25727.94369,-767345.8295.78025595.81 

1160  DATA   -9247486589, 6. 97688E+1 1 , -2. 66192E+13, 3. 94078E-H4 

1170  READ  D(+) 

1180  PRINTER  IS  701 

1190  BEEP 

1200  INPUT  "ENTER  MONTH,  DAY,  AND  TIME  (MM:DD:HH:MM:SS) " , TimeS 

1210  OUTPUT  709;"TD";Time$ 

1220  BEEP 

1230  INPUT  "ENTER  INPUT  MODE  ( 0-3054A, 1 -FILE) *, Ira 

1240  IF  Im-1  THEN 

1250  BEEP 

1260  INPUT  "ENTER  NAME  OF  EXISTING  DATA  FILE"  ,01df i le$ 

1270  PRINT  USING  "20X, ""These  Results  Are  From  Data  File:   "" , 1 OA" ;01df i le$ 

1280  PRINT 

1290  PRINT 

1300  ELSE 

1310  BEEP 

1320  INPUT  "ENTER  NAME  OF  NEW  DATA  FILE?" ,Neuf i leS 

1330  PRINT  USING  "20X,""Data  For  This  Run  Is  Stored  In  Data  File:   "" ,1  Oft" :Newf 

ileS 

1340  PRINT 

1350  PRINT 

1360  END  IF 

1370! 

1380  !   ENTER  THE  CHANNEL  CONFIGURATION. 

1390! 

1400  PRINT  USING  "9X,""THE  FOLLOWING  CHANNEL  CONFIGURATION  WAS  SELECTED: """ 

1410  PRINT 

1420  BEEP 

1430  INPUT  "SELECT  HEATER  TYPE  (O-STRAIGHT.  1 -CURVED) ". I  type 

1440  IF  Itype-0  THEN  PRINT  USING  "12X, " "Heat ing  STRAIGHT  Test  Section. "*" 

1450  IF  I  type- 1  THEN  PRINT  USING  " 12X, ""Heatinq  CURVED  Test  Section. """ 

1460  BEEP 

1470  INPUT  "SELECT  HEATER  POSITION  (0-OUTER,  1 -INNER) ". Ipos 

1480  IF  Ipos-0  THEN  PRINT  USING  " 1 2X , " "Heat ing  OUTER  Plate.*"" 

1490  IF  Ipos-1  THEN  PRINT  USING  " 12X, ""Heating  INNER  Plate. """ 

1500  PRINT 

1510! 

1520  !   ENTERS  FILES  WITH  T/C  CALIBRATION  COEFFICIENTS  BASED  ON  HEATER  TYPE. 

1530! 

1540  IF  Itype-0  THEN  ASSIGN  SFilec  TO  "STRA_COEF" 

1550  IF  Itype-1  THEN  ASSIGN  GFilec  TO  "CURV_C0EF" 

1560  FOR  1-0  TO  76 

1570  ENTER  ©Fi lee: Aa( I , 0 ) , Aa( I . 1 ) , Aa ( I .2) 

1580  NEXT  I 
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1590  ASSIGN  ©Filec  TO  * 

1600! 

1610  !   INSERT  DATA  STORAGE  DISK. 

1620! 

1630  BEEP 

1640  INPUT  "CHANGE  DISK  AND  HIT  ENTER*, Ok 

16S0  IF  Im-1  THEN  ASSIGN  ©File  TO  OldfileS 

1660  IF  Im-0  THEN 

1670  CREATE  BDAT  Newfile$,40 

1680  ASSIGN  ©File  TO  NeufileS 

1690  END  IF 

1700  ENTER  709;Time$ 

1710  CLEAR  709 

1720  PRINT  USING  n20X, "  "Honth,  Day,  and  Time:  ""  .  14ft"  ;Tin»e$ 

1730  PRINT 

1740! 

1750  !   FOR  DATA  RUN  READ  AND  STORE  ALL  RAW  EMF  VALUES. 

1760! 

1770  IF  Int-0  THEN 

1780  OUTPUT  709; "AR  AFOO  AL79" 

1790  OUTPUT  722:*F1  R1  T1  Z1  FL1" 

1800  FOR  1-0  TO  79 

1810  OUTPUT  709;"AS" 

1820  ENTER  722:Emf(I) 

1830  NEXT  I 

1840  OUTPUT  ©File;Emf(*) 

1850  ELSE 

1860  ENTER  ©File;Emf(*) 

1870  END  IF 

1880  OUTPUT  709; "TD" 

1890! 

1900  !   CONVERT  AND  CALIBRATE  ALL  T/C  READINGS  TO  DEGREES  CELSIUS. 

1910! 

1920  PRINT  USING  "9X,""THE  FOLLOWING  DATA  WERE  RECORDED: """ 

1930  PRINT 

1940  FOR  1-0  TO  79 

1950  IF  1-38  OR  1=39  OR  1-77  OR  1-78  THEN 

1960  IF  (Ipos-0  AND  1=38)  OR  (Ipos-1  AND  1-77)  THEN 

1970  Iadd-0 

1980  IF  Ipos-1  THEN  Iadd-39 

1990  PRINT  USING  " 12X, ""Heater  Voltage  (Vh)      =  "",2D.3D,""  ( V) " "" ;Emf (38+Ia 

dd) 

2000  PRINT  USING  " 12X, ""Resistor  Voltage  (Vr)    =  "",2D.3D,nn  ( V) " " " ;Emf (39+Ia 

dd) 

2010  IF  Ipos-0  THEN  PRINT  USING  " 12X, " "Precision  Resistor  (Rpr)  -   2.01729  (Ohm 

s)"u" 

2020  IF  Ipos-1  THEN  PRINT  USING  " 12X, " "Precision  Resistor  (Rpr)  -   2.00839  (Ohm 

s ) "  " " 

2030  PRINT 

2040  END  IF 

2050  ELSE 

2060! 

2070  !   IF  T/C  OUT  OF  RANGE  PRINT  A  QUE. 

2080! 

2090  Aemf-Emf(I) 

2100  IF  Aemf<3.E-4  OR  Aemf>3.E-3  THEN 

2110  T(I)=999.99 

2120  ELSE 

2130  CALL  Tvsv(Emf ( I ) ,Tt ) 

2140  IF  1=79  THEN 
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2150 

2160 

2170 

2180 

2190 

2200 

2210 

2220! 

2230 

2240! 

2250 

2260 
G... 

2270 

(5) 

2280 

1  O  n  n  u 

2290 

T(11) 

2300 

18""" 

2310 

16), T 

2320 

24""" 

2330 

22),T 

2340 

30""" 

2350 

28),T 

2360 

2370! 

2380 

2390! 

2400 

2410 

2420 

2430 

2440 

2450 

2460! 

2470 

ATER. 

2480! 

2490 

2500 

2510 

2520 

2530 

2540 

2550 

2560 

2570 

ii  ■ 

2580 
2590 
2600 
2610 
2620 


T(I)«Tt 

GOTO  2210 

END  IF 

T(I)-FNTcorr(Tt,D 

END  IF 

END  IF 

NEXT  I 

!   PRINT  ALL  T/C  READINGS. 

PRINT  USING  "12X, ""OUTER  PLATE: 
PRINT  USING  "14X,""T/C  Number: 


PRINT  USING 

PRINT  USING 

PRINT  USING 

PRINT  USING 

PRINT  USING 
(17) 
PRINT  USING 

PRINT  USING 
(23) 
PRINT  USING 

PRINT  USING 
(29) 
PRINT 


14X,""Temp  (C): 
'MX, ""T/C  Number: 
14X,""Temp  (C): 
14X,*"T/C  Number: 
14X,""Temp  (C): 
14X,""T/C  Number: 
'14X,""Temp  (C): 
'14X, ""T/C  Number: 
'14X,""Temp  (C): 


12       3       4       5 
,6(3D.DD,2X)":T(0),T(1),T(2),T(3).T(4),T 
7       8       9      10      11 
,6(3D.DD,2X)";T(6),T(7),T(8),T(9),T(10), 
13      14      15      16      17 
,6(3D.DD.2X)*:T(12),T(13),T(14),T(15).T( 
19      20      21      22      23 

'  ,6(3D.DD,2X)";T(18),T(19) ,T(20).T(21 ),T( 
25      26      27      28      29 

' , 6 ( 3D . DD , 2X ) " : T ( 24 ) , T ( 25 ) . T ( 26 ) . T ( 27 ) , T ( 


!   CALCULATE  AVERAGE  TEMPERATURE  OF  OUTER  HEATED  WALL. 

Sum»0. 

IF  Itype-0  THEN 

FOR  1-0  TO  29 

Sum-Sum+T ( I ) 

NEXT  I 

Two-Sum/30 

!   CORRECTION  FOR  BROKEN  T/C  CHANNELS  5,  1G,  21,  AND  26  IN  OUTER  CURVED  HE 


ELSE 

FOR  1-0  TO  29 

IF  1-5  OR  1-16  OR  1-21  OR  1-26  THEN  2530 

Sum-Sum+T ( I ) 

NEXT  I 

Two-Sum/26 

END  IF 

PRINT  USING  "12X, ""OUTER  INSULATION:""" 

PRINT  USING  "14X,"UT/C  Number:     31  (Tinsol) 


32    (Tinso2)        33    (Tinso3) 


PRINT   USING 
Tinso1-T(30) 
Tinso2=T(31 ) 
Tinso3-T(32) 
PRINT 


14X,""Temp    (C) 


,  3  ( 3D .  DD ,  8X  ) " :  T  ( 30  )  ,  T  ( 3 1  )  .  T  ( 32 ) 
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2G30  PRINT  USING  " 12X, " "ORIFICE  TEMP  (Torf)  -  "" ,3D.DD, "  (C)""";T(33) 

2640  PRINT 

2650  PRINT  USING  "12X, ""INLET  TEMPERATURE:""" 

2660  PRINT  USING  "14X,""T/C  Number:    35      36      37      38""" 

2670  PRINT  USING  "14X,""Temp  (C):     "" ,4(3D.DD,2X) " ;T(34) , T(35) ,T(36) ,T(37) 

2680  PRINT 

2690! 

2700  !   CALCULATE  AVERAGE  INLET  TEMPERATURE. 

2710! 

2720  Sum-0. 

2730  FOR  1-34  TO  37 

2740  Sum-Sum+T(I) 

2750  NEXT  I 

2760  Tin-Sum/4 

2770  PRINT  USING  "12X, ""INNER  PLATE:""" 

2780  PRINT  USING  "14X,""T/C  Number:    41      42      43      44      45 

46""" 

2790  PRINT  USING  "14X, ""Temp  (C):     "" ,6(3D.DD,2X) " ; T(40) , T(41 ) ,T(42) ,T(43) ,T( 

44),T(45) 

2800  PRINT  USING  *14X,""T/C  Number:    47      48      49      50      51 

52""" 

2810  PRINT  USING  "14X,""TemP  (C):     " " ,6(3D.DD,2X) " ; T(46) , T(47) , T(48) ,T(49) , T( 

50),T(51) 

2820  PRINT  USING  "14X,""T/C  Number:    53      54      55      56      57 

r  p  «  ■  •! 

2830   PRINT  USING  "14X,""TemP  (C):     "" ,6(3D.DD,2X) " ;T(52) , T(53) , T(54) ,T(55) ,T( 

S6),T(57) 

2840  PRINT  USING  "14X,"*T/C  Number:    59      60      61      62      63 

64""" 

2850   PRINT  USING  "14X,""Temp  (C):     " " ,6(3D.DD,2X) " : T(58) , T(59) , T(60 ) ,T(61 ) ,T( 

62),T(63) 

2860   PRINT  USING  "14X,"*T/C  Number:    65      66      67      68      69  ' 

70""" 

2870  PRINT  USING  "14X,""Temp  (C):    " " ,6(3D.DD,2X) " ; T(64) , T(65) , T(66) ,T(67) ,T( 

68),T(69) 

2880  PRINT 

2890! 

2900   !   CALCULATE  AVERAGE  TEMPERATURE  OF  INNER  HEATED  WALL. 

2910! 

2920   Sum-0. 

2930   IF  Itype-0  THEN 

2940   FOR  1-40  TO  69 

2950  Sum-Sum+T(I) 

2960  NEXT  I 

2970   Tui-Sum/30 

2980! 

2990   !   CORRECTION  FOR  BROKEN  T/C  CHANNEL  47  IN  INNER  CURVED  HEATER. 

3000! 

3010  ELSE 

3020  FOR  1-40  TO  69 

3030  IF  1-47  THEN  3050 

3040  Sum-Sum+T(I) 

3050  NEXT  I 

3060  Tun-Sum/29 

3070  END  IF 

3080   PRINT  USING  " 1 2X, " "INNER  INSULATION:""* 

3090  PRINT  USING  "14X,""T/C  Number:     71  (TinsiD    72  (Tinsi2)    73  (Tinsi3)" 

ii  n 

3100   PRINT  USING  " 1 4X. " "Temp  (C):        " " . 3(3D . DD , 8X ) " ; T ( 70 ) , T ( 7 1  )  , T ( 72 ) 
3110   Tinsi1-T(70) 
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3120 

3130 

3140 

3150 

3160 

3170 

3180 

3190! 

3200 

3210! 

3220 

3230 

3240 

3250 

32G0 

3270! 

3280 

3290! 

3300 

3310! 

3320 

3330! 

3340 

3350! 

3360 

3370! 

3380 

3390 

3400 

DD,"" 

3410 

DD , "  " 

3420 

DD,"" 

3430 

DD,"" 

3440 

DD,"" 

3450 

DD,"" 

3460! 

3470 

3480! 

3490 

3500 

3510 

3520 

3530 

DD,"" 

3540! 

3550 

3560! 

3570 

3580 

DD,"" 

3590! 

3G00 

3510! 

3G20 

3G30 


Tinsi2«T(71) 

Tinsi3-T(72) 

PRINT 

PRINT  USING  "12X, ""OUTLET  TEMPERATURE: " "" 

PRINT  USING  "14X.""T/C  Number:    74      75      76      77""" 

PRINT  USING  "14X,""TemP  (C):     "" , 4 (3D . DD,2X) * ; T( 73) , T ( 74) , T ( 75) , T ( 76 ) 

PRINT  USING  "«,*" 

!   CALCULATE  AVERAGE  OUTLET  TEMPERATURE. 

Sum«0. 

FOR  1-73  TO  76 

Sum-Sum+T( I ) 

NEXT  I 

Tout-Sum/4 

!   CALCULATE  DIFFERENCE  BETWEEN  CHANNEL  INLET  AND  OUTLET  TEMPERATURES. 

Tdiff-Tout-Tin 

!   CALCULATE  CHANNEL  BULK  TEMPERATURE. 

Tblk-(Tin+Tout)/2 

!   PRINT  CALCULATED  TEMPERATURES. 

PRINT  USING  "9X,""THE  FOLLOWING  TEMPERATURES  WERE  CALCULATED: " "" 

PRINT 

PRINT  USING  "12X, ""Average  Outer  Wall  Temperature  (Two)  ■  "".3D. 

(C)""";Two 

PRINT  USING  "12X, ""Average  Inner  Wall  Temperature  (Tun)  -  "".3D. 

(C)""";Tui 

PRINT  USING  "12X, ""Average  Outlet  Temperature  (Tout)  -  "".3D. 

(C)""";Tout 

PRINT  USING  "12X, ""Average  Inlet  Temperature  (Tin)  -  ".3D. 

(C)""";Tin 

PRINT  USING  "12X, ""Channel  Inlet  and  Outlet  Temp  Diff  (Tdiff)     -  ",3D. 

(C)""";Tdiff 

PRINT  USING  "12X, ""Average  Bulk  Temperature  (Tbllc)  -  ",3D. 

(C)""";Tblk 

!   CALCULATE  DELTA  T  (Tdel)  FOR  USE  IN  CALCULATING  H  LATER. 

Tdelo-Tuo-Tblk 

Tdeli-Twi-Tblk 

IF  Ipos-0  THEN  Tdel-Tdelo 

IF  Ipos-1  THEN  Tdel-Tdeli 

PRINT  USING  "12X,""Mean  Temperature  Difference  (Tdel)  ■  "".3D. 

(C)""";Tdel 

!   CALCULATE  DIFFERENCE  BETWEEN  OUTER  AND  INNER  WALL. 

Twdiff-ABS(Tuo-Tuji) 

PRINT  USING  "12X, ""Outer  and  Inner  Wall  Temp  Difference  (Tudiff)   -  "".3D, 

(C)""";Tuidiff 

!   CALCULATE  LOCAL  CHANNEL  TEMPERATURES  AT  INLET.  MIDDLE  AND  OUTLET. 

Tcin=Tin+( .5*Tdiff/12) 
Tcmid=Tin+(4.S*Tdiff/12) 
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3640      Tcout-Tin+U1.5*Tdiff/12> 

3650   PRINT  USING  * 12X, ""Average  Local  Channel  Temp  at  Inlet  (Tcin)     -  "",3D. 

DD,""  (C)""";Tcin 

3660   PRINT  USING  " 12X, " "Average  Local  Channel  Temp  at  Middle  (Tcmid)    -  "",3D. 

DD,""  (C)""";Tcmid 

3670   PRINT  USING  " 12X, ""Average  Local  Channel  Temp  at  Outlet  (Tcout)    -  "",3D. 

DD,""  (C)""";Tcout 

3680! 

3690   !   CALCULATE  AVERAGE  LOCAL  HEATED  UALL  TEMPERATURE  AT  INLET. 

3700! 

3710   Sum-0. 

3720   Iadd-0 

3730   IF  Ipos-1  THEN  Iadd-40 

3740! 

37S0   !   CORRECTION  FOR  BROKEN  T/C  CHANNELS  B,  AND  47  CURVED  HEATER  ONLY. 

3760! 

3770   IF  Itype-1  THEN 

3780   FOR  1-0  TO  6 

3790   IF  (Ipos-0  AND  1-5)  OR  (Ipos-1  AND  1-6)  THEN  3810 

3800   Sum-Sum+T(I+Iadd) 

3810   NEXT  I 

3820   Twin-Sum/6 

3830  ELSE 

3840   FOR  1-0  TO  7 

3850   IF  (Ipos-0  AND  1-7)  OR  (Ipos-1  AND  1-6)  THEN  3870 

3860   Sum-Sum+T(I+Iadd) 

3870   NEXT  I 

3880   Tuin-Sum/7 

3890  END  IF 

3900   PRINT  USING  "12X,""Avg  Local  Heated  Wall  Temp  at  Inlet  (Twin)     -  "".3D. 

DD,""  (C) ;Twin 

3910! 

3920   !   CALCULATE  AVERAGE  LOCAL  HEATED  HALL  TEMPERATURE  NEAR  MIDDLE. 

3930! 

3940   Sum-0. 

3950! 

3960   !   CORRECTION  FOR  BROKEN  T/C  CHANNEL  16  CURVED  HEATER  ONLY. 

3970! 

3980   IF  Itype-1  AND  Ipos-0  THEN 

3990  FOR  1-9  TO  15 

-4000   IF  1-10  THEN  4020 

4010   Sum-Sum+T(I) 

4020   NEXT  I 

4030   Tumid=Sum/6 

4040   ELSE 

4050   FOR  1-9  TO  16 

4060   IF  (Ipos-0  AND  1-10)  OR  (Ipos-1  AND  1-9)  THEN  4080 

4070   Sum-Sum+T(I+Iadd) 

4080   NEXT  I 

4090   Twmid-Sum/7 

4100   END  IF 

4110   PRINT  USING  "12X,""Avg  Local  Heated  Wall  Temp  Near  Middle  (Tumid)  -  "".3D, 

DD,""  (C)""";Tu)mid 

4120! 

4130   !   CALCULATE  AVERAGE  LOCAL  HEATED  WALL  TEMPERATURE  AT  OUTLET. 

4140! 

4150   Sum-0. 

4160! 

4170   !   CORRECTION  FOR  BROKEN  T/C  CHANNEL  26  CURVED  HEATER  ONLY. 

4180! 
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4190 

4200 

4210 

4220 

4230 

4240 

4250 

4260 

4270 

4280 

4290 

4300 

4310 

DD,"" 

4320 

4330 

4340! 

4350 

43G0! 

4370 

4380 

4390 

4400 

4410 

4420 

4430 

4440! 

4450 

44G0! 

4470 

4480 

4490 

4500 

4510! 

4520 

4530! 

4540 

4550 

4550 

4570 

4580 

4590 

4600 

4610! 

4620 

4630! 

4640 

4650 

4660 

4670 

4680 

4690 

4700 

4710! 

4720 

4730! 

4740 

4750 

4760! 

4770 


IF  Itype-1  AND  Ipos-0  THEN 

FOR  1-23  TO  29 

IF  1-26  THEN  4230 

Sum-Sum+T ( I ) 

NEXT  I 

Twout-Sum/6 

ELSE 

FOR  1-23  TO  29 

Sum-Sum+T(I+Iadd) 

NEXT  I 

Twout-Sum/7 

END  IF 

PRINT  USING  "12X,"Avg  Local  Heated  Wall  Temp  at  Outlet  (Tuout)    -  ".3D. 

(C)"";T«uout 

PRINT 

PRINT 

!   ENTER  ORIFICE  PRESSURE  DATA. 

BEEP 

INPUT  "ENTER  PATH  (inHg),  DPMQnH20).  P1M(inH20).  RE  1  "  .  Patm.  Dpra.P1  m.Rel 

PRINT  USING  "9X,""THE  FOLLOWING  ORIFICE  DATA  WERE  ENTERED:""" 

PRINT 

PRINT  USING  "T2X,""Patm(inHg)      DPM(inH20)     P1MUnH20)      RET" 

PRINT  USING  "14X,2(2D.DD,10X),(2D.DD,7X),5D.D*;Patm.Dpm,P1m,Re1 

PRINT 

!   CONVERT  PRESSURE  READINGS  TO  SI  UNITS  AND  P1M  TO  ABSOLUTE. 

Patm-Patm+3376.8 
Pdel-Dpm*248.84 
P1-Patm-(P1m*248.84) 
Torf-T(33) 

!   PHYSICAL  PROPERTIES  AND  CONSTANTS. 

R-286  987  'FOR  AIR 

Rho-Pl'/(R*(Torf+273.15))  (DENSITY  OF  AIR  AT  ORIF  PLATE. 

Gamma-1 .40  !FOR  AIR. 

Cp-1006  !Cp  GOOD  FOR  Tblk  BETWEEN  12  AND  33  (DEG  C). 

Gc-1.0 

Siqma-5.669E-8 

Ecu-. 12   !EMISSIVITY  OF  COPPER. 

!   DIMENSIONS  OF  CHANNEL. 

Dc-. 006350 

Wid-.2540 

Ri-.297 

Puet-2*(Dc+Wid) 

Ac-Dc+Wid 

Dhd-4*Ac/Puet 

Apl-. 0774192 

!   PROPERTIES  OF  INSULATION. 

Xins-.0127 
Kins-3.8E-2 

!   ENTER  ORIFICE  CONFIGURATION. 


87 


4780! 

4790   BEEP 

4800   INPUT  "SELECT  DIAMETER  OF  ORIFICE  (0-.5334,  1-1.0755  ( inches) )" .Size 

4810   IF  Size-0  THEN  Dorf-. 013548 

4820   IF  Size-1  THEN  Dorf-. 027318 

4830! 

4840   !   DIMENSIONS  OF  PIPE  AND  ORIFICE  PLATE. 

4850! 

4860   A-(PI*Dorf *2)/4 

4870   Dpipe-. 051772 

4880  Apipe-PI*Dpipe"2/4 

4890   Beta-Dorf/Dpipe 

4900   PRINT  USING  " 12X, ""Dorf (m)        A(m*2)         Beta""" 

4910  PRINT  USING  "12X, (Z.6D.6X) , (Z.3DE.6X) ,Z.4D" ;Dorf , A, Beta 

4920  PRINT 

4930  PRINT 

4940! 

4950   !   CORRELATION  FOR  EXPANSION  FACTOR  BASED  ON  ID  AND  1 /2D  TAPS. 

49G0! 

4970   Y-1-(.41  +  .35*Beta'4)MPdel/(Gamroa+P1)) 

4980! 

4990   !   CALCULATION  OF  FLOW  COEFFICIENT  (K)  IN  SI  UNITS. 

5000! 

5010  B-. 0002+2. 794E-5/Dpipe+( .0038+1 . 016E-S/Dpipe)*(Beta*2+( 16.5+1 .968504E+2*Dp 

ipe)*Beta*16) 

5020! 

5030   !   K1  AND  K2  ARE  DUMMY  VARIABLES  TO  CALCULATE  Ko. 

5040! 

5050   K1-.6014-5.3974064E-3/Dpipe'.25 

5060  K2-(.376+2.8971138E-2/Dpipe'.25)*(1.6129E-7/(Dpipe-2*Beta,2+6.3SE-5*Dpipe) 

+Beta*4+1 .5*Beta*16) 

5070   Ko«K1+K2 

5080  K-Ko+1000*B/Re1'.S 

5090! 

5100   !   Mu  CORRELATION  GOOD  FOR  Torf  BETWEEN  17  AND  41  DEGREES  CELSIUS. 

5110! 

5120  Mu-4.6971429E-8*Torf+1 .7194722E-5 

5130! 

5140   !   CALCULATE  Mdot  AND  REpipe  AND  COMPARE  TO  PREDICTED  RE1 . 

5150! 

5160  Mdot-K*A*Y*(2*Gc+Rho*Pdel)* .5 

5170  Repipe-(Mdot*Dpipe)/(Mu*Apipe) 

5180  Diff-(ABS(Re1-Repipe)/Re1 )*100 

5190   Re1-Repipe 

5200   IF  DifO.001  THEN  5080 

5210   PRINT  USING  "9X,""THE  FOLLOWING  DATA  WERE  CALCULATED:""" 

5220  PRINT 

5230   PRINT  USING  " 12X, ""Or  if  ice  Expansion  Factor  (Y)    -   "",Z.4D,";Y 

5240  PRINT  USING  "12X, ""Or  if  ice  Flow  Coefficient  (K)    -   "'\Z.4D,":K 

5250   PRINT  USING  "  12X,  ""Density  Based  on  Torf  (Rho)     -   "".Z.4D,""     (Icg/m* 

3)""" ;Rho 

5260   PRINT  USING  " 12X, " "Viscosi ty  Based  on  Torf  (Mu)    -  "".2Z.3DE,""  (kg/m.s) 

;Mu 

5270! 

5280   !   Mu  CORRELATION  GOOD  FOR  Tblk  BETWEEN  17  AND  41  DEGREES  CELSIUS. 

5290! 

5300      Mu-4.6971429E-8*Tblk+1 .7194722E-5 

5310      Red=(Mdot*Dc)/(Mu*Ac) 

5320      Rehd-(Mdot*Dhd)/(Mu+Ac) 

5330! 
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5340 

5350! 

5360 

5370 

5380 

5390 

5400 

5410! 

5420 

5430! 

5440 

5450! 

54G0 

5470! 

5480 

5490 

5500 

5510 

5520 

5530 

5540 

5550 

5560 

5570! 

5580 

5590! 

5600 

5610! 

5620 

5630! 

5640 

5650 

5660! 

5670 

5680! 

5690 

5700 

5710 

5720 

5730 

5740 

5750 

5760 

5770 

5780! 

5790 

5800! 

5810 

5820 

"■•;M 

5830 

" " ;  Ka 

5840 

5850 

5860 

5870 

5880 

5890 

5900 

5910 


!   CALCULATION  OF  POWER  INTO  HEATER  PLATE. 

Iadd-0 

IF  Ipos-1  THEN  Iadd-39 

IF  Ipos-0  THEN  Rpr-2. 01729 

IF  Ipos-1  THEN  Rpr-2. 00839 

Qp-(En>f(38+Iadd)*Emf  (39+Iadd) ) /Rpr 

!   CALCULATION  OF  HEAT  INPUT  TO  AIR. 

Oair-Mdot+Cp* (Tout-Tin) 

!   CALCULATE  LOCAL  HEAT  TRANSFER  COEF.  AND  NUSSELT  NUMBER. 

Kai r-7. 7257 143E-5*Tcin+. 024 165836 

Hin-Qair/(Apl*(Twin-Tcin) ) 

Nuin-(Hin*Dhd)/Kair 

Kair-7 . 72571 43E-5*Tcm id+ . 0241 65836 

Hmid-Qair/(Apl*(TuKnid-Tcmid)  ) 

Numid-(Hmid*Dhd)/Kair 

Kai r-7. 7257 143E-5*Tcout+. 024 165836 

Hout-Qair/(Apl*(Twout-Tcout) ) 

Nuout-(Hout*Dhd)/Kair 

!   Kair  CORRELATION  GOOD  FOR  Tblk  BETWEEN  17  AND  41  DEGREES  CELSIUS. 

Kai r-7. 7257 143E-5*Tb Ik +.024 165836 

!   CALCULATE  AVERAGE  HEAT  TRANSFER  COEF.  AND  NUSSELT  NUMBER. 

Havg-Qair/(APl*Tdel) 
Nuavg- ( Havg*Dhd ) /Ka  i  r 

!   CALCULATE  HEAT  LOSSES. 

Qlo1-(Kins*Apl*(Tinso1-Tinso2) )/Xins 
Qlo2-(Kins*Apl«-(Tinso2-Tinso3)  )/Xins 
Qlo-(Qlo1+Qlo2)/2 

Qli1-(Kins*Apl*(Tinsi1-Tinsi2) )/Xins 
Qli2-(Kins*Apl*(Tinsi2-Tinsi3) )/Xins 
Qli-(Qli1+Qli2)/2 

Rr-2.0*((1-Ecu)/(Apl*Ecu))+(1 .O/Apl) 
Qr-(Sigma*(ABS((Two+273.1S)"4-(Tui+273.15)'4)) )/Rr 
Odel-Qp-Qair-Qlo1-Qr-Qli1 

!   CALCULATE  DEAN  NUMBER. 

De-.5*Red*((Dc/2)/Ri)'.5 

PRINT  USING  "12X, "'Viscosity  Based  on  Tblk  (Mu)  -  "".2Z.3DE, 
u 

PRINT  USING  "12X, ""Therm  Cond  Based  on  Tblk  (Kair)  -  "",Z.3DE, 
ir 

PRINT 

PRINT  USING  "12X,""Patm   -  "".6D.DD,"" 

PRINT  USING  "12X,""Pdel   -    "",4D.DD, 

PRINT  USING  "12X,""P1     -  "".6D.DD,"" 

PRINT  USING  "12X,""Mdot   -    "",Z.4D, 

PRINT  USING  M12X,""Repipe  -  "",6D 

PRINT  USING  "12X,""Red    -  "" ,6D 

PRINT  USING  "12X.""Rehd   -  "".6D 


"  (kg/ro.s) 
"  (W/m.K)" 


D" : Repipe 
,D":Red 
,D":Rehd 


(N/ro'2)"" 
"  (N/m"2) 
(N/tn*2)"" 
(kg/s)" 


:Patm 
"":Pdel 

;P1 

"  :Mdot 
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5920! 

5930 

5940! 

5950   IF  Itype-1  THEN  PRINT  USING  "12X, 

5960  PRINT 

5970  PRINT  USING  "12X,""Qp 

atts)"";Qp,Qair 

5980  PRINT  USING  " 12X, " "Qlo1 

<Watts)"";Qlo1,Qli1 

5990  PRINT  USING  "12X,""Qlo2 

(Watts)"";Qlo2,Qli2 

G000  PRINT  USING  "12X,""Qlo 

(Watts)"";Qlo.Qli 

6010  PRINT  USING  "12X,""Qr 

(Watts)"";Qr,Qdel 

6020   PRINT 

6030   PRINT  USING  "12X,"Hin 

6040   PRINT  USING  "12X,"Hmid 
mid.Numid 

6050   PRINT  USING  "12X,"Hout 
out , Nuou t 

6060   PRINT  USING  "12X,""Havg 
avg,Nuavg 
607Q   PRINT 
PRINT 


ONLY  PRINT  DEAN  NUMBER  WHEN  HEATING  CURVED  SECTION. 

De     -   ",4D.D";De 
",3D.3D,"  (Watts)", 7X,"Qair  - 
",2D.3D,"  (Watts)"", 7X,""Qli1  ■ 
"",2D.3D,""  (Watts)"", 7X,""Qli2  ■ 
"",2D.3D,""  (Watts)"", 7X,""Qli   ■ 
'".2D. 3D,"  (Watts)"", 7X,*"Qdel  ■ 


.3D.DD," 
,3D.DD,"' 
,3D.DD," 


(U/m*2C)' 

(W/m'2C) 

(W/m*2C) 


,3D.DD,""  (W/m*2C) 


,7X,""Nuin 
,7X,"Numid 
,7X,""Nuout 
,7X,""Nuavg 


,3D.3D,""  (W 

",2D.3D," 

"",2D.3D," 

",2D.3D," 

"",2D.3D," 

"",3D.DD,";H 
"",3D.DD,";H 
"",3D.DD,";H 
",3D.DD,";H 


6080 

6090 

6100 

6110 

6120 

6130 

6140 

6150 

6160 

6170 

6180 

6190 

6200 

6210! 

6220 

6230! 

6240 

6250 

6260 

6270 

6280 

6290 

6300 

6310 

6320! 

6330 


BEEP 

INPUT  "WILL  THERE  BE  ANOTHER  RUN?  ( 1 -Y ,0-N) " ,Go_on 

IF  Go.onOO  THEN 

PRINT  USING  "$,*" 

IF  Im-0  THEN  WAIT  600 

GOTO  1700 

END  IF 

PRINT  USING  "8X,""END  OF  RUN"" 

PRINT  USING  "«>,#" 

PRINTER  IS  1 

ASSIGN  ©File  TO  * 

END 

!   CONVERTS  EMF  TO  DEGREES  CELSIUS. 

SUB  Tvsv/(V,T) 

COM  /Co/  D(7),Aa(76,2) 

Sum"0 

FOR  1-0  TO  7 

Sum=Sum+D(I)*V*I 

NEXT  I 

T  =  Sum 

SUBEND 

!   CALIBRATES  T/C  READINGS. 


6340! 

6350  DEF  FNTcorr(T,I) 

6360  COM  /Co/  D(7),Aa(76,2) 

6370  Tc-Aa(I.O) 

6380  FOR  J-1  TO  2 

6390  Tc-Tc+Aa(I,J)*T*J 

6400  NEXT  J 

6410  RETURN  Tc 

6420  FNEND 
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APPENDIX  B:   SAMPLE  CALCULATIONS 

Figure  36  below,  shows  the  major  heat  transfer 
components  for  each  of  the  test  sections.   The  sample 
calculations  that  follow,  demonstrate  the  methods  used  by 
the  computer  to  calculate  these  components  as  well  as  the 
Reynolds  number,  average  heat  transfer  coefficient,  and 
average  Nusselt  number  for  each  set  of  data.   Part  A  of 
Appendix  B  contains  a  sample  printout  of  a  data  run, 
followed  by  a  summation  of  the  data.   The  sample 
calculations  are  for  the  concave  curved  section,  but  those 
for  the  straight  section  are  similar. 


Qp 


HEATED  COPPER  WALL 


^  Qlo 


~1 


I 


Qr 


l 


Q  a  ir, 


L 


7TrT?ri 


1 


Figure  36.   Energy  Balance  in  Test  Section 
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1.   SAMPLE  CALCULATION  DATA 

THE  FOLLOWING  CHANNEL  CONFIGURATION  WAS  SELECTED: 

Heating  CURVED  Test  Section. 
Heating  OUTER  Plate. 

Month,  Day,  and  Time:  10:17:13:47:11 

THE  FOLLOWING  DATA  WERE  RECORDED: 


Heater  Voltage 

(Vh) 

-  61  . 

337  (V) 

Resistor  Voltag 

e  (Vr) 

«  5. 

195  (V) 

Precision  Resis 

tor  (Rpr 

)  -  2. 

01729  (Ohms) 

OUTER  PLATE: 

T/C  Number: 

1 

2 

3 

4 

5 

6 

Temp  ( C )  : 

42.75 

41  .57 

44.39 

41  .23 

44.55 

40.59 

T/C  Number: 

7 

8 

9 

10 

11 

12 

Temp  ( C ) : 

46.54 

42.61 

43.51 

43.47 

44.77 

43.86 

T/C  Number: 

13 

14 

15 

16 

17 

18 

Temp  (C) : 

43.40 

44.18 

44.10 

41.97 

38.06 

45.73 

T/C  Number: 

19 

20 

21 

22 

23 

24 

Temp  ( C) : 

45.79 

48.45 

49.18 

33.83 

45.32 

45.51 

T/C  Number: 

25 

26 

27 

28 

29 

30 

Temp  ( C)  : 

45.30 

46.47 

30.31 

45.80 

45.98 

47.28 

OUTER  INSULATION: 

T/C  Number: 

31  (Ti 

msol  ) 

32  (Tinso2) 

( Tinso 

3) 

Temp  (C) : 

49 

.58 

37.29 

29.36 

ORIFICE  TEMP  (T 

orf)  - 

24.16  ( 

C) 

INLET  TEMPERATURE: 

T/C  Number: 

35 

36 

37 

38 

Temp  ( C) : 

20.55 

20.48 

20.47 

20.43 

INNER  PLATE: 

T/C  Number: 

41 

42 

43 

44 

45 

46 

Temp  (C) : 

21  .43 

21  .31 

21  .16 

21  .14 

21  .05 

21.92 

T/C  Number: 

47 

48 

49 

50 

51 

52 

Temp  ( C ) : 

21  .03 

22.39 

21  .10 

21.05 

21  .38 

21  .46 

T/C  Number: 

53 

54 

55 

56 

57 

58 

Temp  ( C) : 

21  .23 

21.27 

21  .17 

21.19 

21  .15 

21  .35 

T/C  Number! 

59 

60 

61 

62 

63 

64 

Temp  (C) : 

21.25 

21.56 

21  .47 

21  .74 

21  .66 

22.18 

T/C  Number: 

65 

66 

67 

68 

69 

70 

Temp  (C): 

22.05 

21.91 

21.85 

21  .78 

21.74 

21  .73 

INNER  INSULATION: 

T/C  Number: 

71  (T 

insi 1 ) 

72  (Tinsii 

( Tinsi 

3) 

Temp  (C) : 

21 

.53 

21  .49 

21  .32 

OUTLET  TEMPERATURE: 

T/C  Number: 

74 

75 

76 

77 

Temp  ( C) : 

24.29 

24.38 

24.38 

24.25 
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THE  FOLLOWING  TEMPERATURES  WERE  CALCULATED 


Average  0 
Average  I 
Average  0 
Average  I 
Channel  I 
Average  B 
Mean  Temp 
Outer  and 
Average  L 
Average  L 
Average  L 
Avg  Local 
Avg  Local 
Avg  Local 


uter  Wall  Te 
nner  Wa 1 1  Te 
utlet  Temper 
nlet  Tempera 
nlet  and  Out 
ulk  Temperat 
erature  Diff 

Inner  Wall 

ocal  Channel 

ocal  Channel 

ocal  Channel 

Heated  Wall 

Heated  Wall 

Heated  Wall 


mperature  (Two)  =  44.76  (C) 

mperature  (Tui)  ■  21.46  (C) 

ature  (Tout)  =  24.32  (C) 

ture  (Tin)  =  20.48  (C) 

let  Temp  Diff  (Tdiff)  =  3.84  (O 

ure  (Tblk)  =  22.40  <C) 

erence  (Tdel)  ■  22.36  (C) 

Temp  Difference  (Twdiff)  =  23.30  (C) 

Temp  at  Inlet  (Tcin)  =  20.64  (C) 

Temp  at  Middle  (Tcmid)  =  21.92  (C) 

Temp  at  Outlet  (Tcout)  =  24.16  (C) 

Temp  at  Inlet  (Twin)  =  43.51  (C) 

Temp  Near  Middle  (Tumid)  =  43.50  (O 

Temp  at  Outlet  (Tuiout)  =  46.06  (C) 


THE  FOLLOWING  ORIFICE  DATA  WERE  ENTERED: 


Pa  tin  (  inHg) 
29.79 


DPM(inH20) 
18.58 


P1M(  inH20) 
12.55 


RE1 
15500 


Dorf (m) 
0.027318 


A(m*2) 
5.861E-04 


Beta 
0.5277 


THE  FOLLOWING  DATA  WERE  CALCULATED: 

Orifice  Expansion  Factor  (Y) 
Orifice  Flow  Coefficient  (K) 
Density  Based  on  Torf  (Rho) 
Viscosity  Based  on  Torf  (Mu) 
Viscosity  Based  on  Tblk  (Mu) 
Therm  Cond  Based  on  Tblk  (Kair) 


0.9852 

0.6323 

1.1424  (kg/m*3) 

18.330E-06  (kg/m.s) 

18.247E-06  (kg/m.s) 

2.590E-02  (W/m.K) 


Patm 

=    100594.87    (N/nT 

'2) 

Pdel 

4623.45    (N/nT 

'2) 

P1 

=      97471  .93    (N/nT 

'2) 

Mdot 

0.0375    (kg/s) 

Rep ipe 

=      50352.2 

Red 

8097.1 

Rehd 

=      15799.1 

De 

418.6 

Qp 

■ 

157.971     (Watts) 

Oair    = 

144.950 

(Watts) 

Qlol 

* 

2.847    (Watts) 

01 1 1    = 

.009 

(Watts) 

Qlo2 

* 

1.838    (Watts) 

Qli2   = 

.038 

(Watts) 

Qlo 

- 

2.342    (Watts) 

Qli      = 

.023 

(Watts) 

Qr 

= 

.751     (Watts) 

Qdel    = 

9.415 

(Watts) 

Hin 

= 

81 .90    (W/m'2C) 

Nu  i  n 

=      39.39 

Hmid 

= 

86.78    (W/m*2C) 

Numid 

=      41.58 

Hout 

= 

85.51     (W/m*2C) 

Nuout 

=      40.70 

Havg 

= 

83.75    (W/m'2C) 

Nuavg 

=      40.07 

93 


VH 
VPR 

RPR 

P  . 
atm 

AP 


air 

K. 
ins 


AX 


ins 


cu 


=  61.337  V 
=  5.195  V 
=  2.0173  Q 
=  29.79  in.  Hg 


=  18.58  in.  H20 


pl 

= 

12.55 

in.  H20 

T. 
in 

= 

20.48 

C 

T 
out 

= 

24.32 

C 

T 
wo 

= 

44.76 

C 

T  . 

Wl 

= 

21.49 

C 

T. 
ins , 

il 

= 

21.53 

C 

T. 
ins, 

i2 

= 

21.49 

C 

T. 
ins , 

01 

= 

49.58 

C 

T- 

ins, 

02 

= 

37.29 

C 

Torf 

= 

24.16 

K  . 

air 

= 

25.90 

x  10"3  W/mC 

cMl, 

= 

1006  J/Kg  K 

=  4.6971  x  10  8  x  T  .   +  1.71947  x  10 

air 

3.80  x  10~2  W/mC 

=  0.0127  m 

=  0.12 

=  5.669  x  10~8  W/m2K4 

=  .5277 

=  1.40 

=  1  Kg  m/Ns2 


-5 
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R 

= 

286.98  Nm/Kg  K 

F 
wo-wi 

= 

1.0 

APL 

= 

.07742  m2 

A 
c 

= 

.0016  m2 

A  . 
pipe 

= 

.00211  m2 

D 
c 

= 

.00635  m 

D  . 
pipe 

= 

.0518  m 

Dorf 

= 

.027318  m 

P 
wet 

= 

.5207  m 

Ri 

= 

.297  m 

2.   TEMPERATURE  CALCULATIONS 

a.   Bulk  Temperature  (T, ,. ) 


TKifc  =  T'"  *  T°"t  ,  20.48  ♦  24.32  . 
blk        2  2 


b.   Mean  Temperature  Difference  (AT) 


AT  =  T    -  T.  ,,  =  44.76  -  22.40  =  22.36 
wo     blk 


c.   Temperature  Difference  Between  Outer  and  Inner  Wall 
(Twdiff} 


Tr,.*<**  =  T    -  T  .  =  44.76  -  21.49  =  23.27  C 
wdiff     wo     wi 
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3.   POWER  CALCULATIONS 


a.   Power  Supplied  (Q_) 


Q   =  !"!™  .  (61.337)(5.195)  .  157.97  „ 
P      RpR         (2.0173) 


b.   Heat  Lost  Through  the  Back  of  the  Outer  Plate  (Qlo) 


Tins,ol    Tins,o2   =  (49.58  -  37.29) 

10    (AXins)/(Kins)  (APL}    (0. 0127) / (3. 80xl0~2)  (.07742) 

=  2.85  W 


c.   Heat  Lost  Through  the  Back  of  the  Inner  Plate  (Q.) 


Tins,il    Tins,i2   _       (21.53  -  21.49) 


11    (AXins)/(Kins)  (APL}    (0. 0127) / (3. 80x10  2)  (.07742) 
=  .01  W 


d.   Heat  Radiated  (Q  ) 

(1)  Radiation  Resistance  (RR) 


1  -  e 
RD  =  2  x  (u -^)  + 


R        APL  ecu    APL  Fwo-wi 


\.ZT~    "     1]      =     X     rAnTrS      ["To    "     H      =     202.36    R|     2 


APL       £cu  °'07742       *12 
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(2)  Heat  Radiated  (Q  ) 


4       4 
a  (T     -  T  .  ) 

°r  =  — ^ ~       T  [=1  * 


(5.669x10  8)  [  (44.76  +  273)4  -  (  21 . 49  +  273) 4] 
202. 36 

=  0.75  W 

4.   MASS  FLOW  RATE  CALCULATIONS 
a.   Pressure  Conversions 

(1)  Patm=  29.79  in  Hg  x  3376.8  =  100594.87  N/m2 

(2)  AP   =  18.58  in  H20  x  248.84  =  4623.5  N/m2 

(3)  P.   =  (12.55  in  H20  x  248.84)  -  100594.87 

=  97471.93  N/m2 


Density  of  Air  (p_,-„) 

air 


Pl  (97471.93) ,     ,.0A    Vri/m3 

air    "    R-TT^   =    (286.  98 )  (24.  16    +      273)'    U14U    Kg/m 


c.       Expansion    Factor     (Y) 


Y    =    1    -    [.41    +    .3564]    ~- 

y    1 


4  4623    5 

=    1    -    [.41    +    .35(.5277)4] 


(1.40)  (97471.93) 
=    .9852 
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d.   Area  of  Orifice  (A) 


A  = 


»  (D  f)2         tt(.027318)2 

orf  ' 


=  .005861  m' 


e.   Mass  Flow  Rate  (m) 


m  =  YKA    2g   p  .   AP 
^c   air 


=  (.9852)K(. 0005861) 


x    2(1)  (1.1424)  (4623.5) 
m  =  .05935  K 

Iterating: 

Assume  a  Pipe  Reynolds  number  =  50000 

Obtain  a  value  for  K,  the  flow  coefficient,  from 

reference  30.  K  =  .6325 
Solve  for  m.  m  =  .0375  Kg/s 
Solve  for  new  Repipe  =   m  Dpipe      (.q375)  (.0518) 


ApipeMair   (• 00211)  (18. 33x10  6) 

Re  .    =  550225 
pipe 

Check  convergence  and  repeat  process  if  necessary. 
(Convergence  if  difference  less  than  .001) 
m  =  .0375  Kg/s 
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5.   REYNOLDS  NUMBER  CALCULATIONS 


a.   Re     =  J^PiPi_  =    (.0375X.0518)      .  50225 
Apipeuair    (-00211)  (18. 33X10"6) 


b. 

Red 

m   D 

c 

Ac    ^air 

r  . 

Rp.     . 

1   Dhd 

(.0375)  (.00635)    =  8Q97>1 
(.0016)  (18.25x10   ) 


hd    A   y  . 
c   air 

^s       n      4  X   Ac    (4)  (.0016)     mooo 
U)   °hd  =   Pwo.    =    .5207     =  -01229  m 

R     .  m   °hd    =   (.0375)  (.01229)     =  1579g>1 

A  p  .      (.0016)  (18. 25x10   ) 
c  ai  r 


6.  HEAT  CONVECTED  TO  AIR  CALCULATION 

a.   Q.*  «  mC  fT.-T,  1  =  (.0375)  (1006)  (24. 32-20. 48)=144. 95  W 

air     p   out   in 

7.  AVERAGE  HEAT  TRANSFER  COEFFICIENT  CALCULATION 


h     °air   .      144.95       .  p-a  79  w/m2r 
h  "  A^Tf  "  (.07742)  (22.36)  "  83*72  W/m  C 


8.   AVERAGE  NUSSELT  NUMBER  CALCULATION 


551  _  =  l^.   =  (83.73)  (.01229)  =  39<73 
hd     air     25.90  x  10~3 
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9.   DEAN  NUMBER  CALCULATION 


Dc/2     8097.1    /. 00635/2   _  ,1Q  c 
TT      =   — \/   .297     "  418*6 


100 


APPENDIX  C:   EXPERIMENTAL  UNCERTAINTY 

The  uncertainties  for  the  major  variable  in  the 
experiments,  were  calculated  in  accordance  with  the  method 
described  by  S.  Kline  and  F.  McClintoch  [Ref.  33].   The 
estimates  of  the  uncertainty  in  the  measured  quantities  were 
made  conservatively.   As  a  result,  there  is  considerable 
confidence  in  the  uncertainties  as  calculated. 

The  following  equations  were  used  to  calculate  the 
uncertainties : 


,       ,D  2     dP,  2     dT   f  2 

1        orf 


1/2 


(2)   ^ 
m 


[<<^>2  +  <f)2  ♦  <^)2) 


+  (1/4  ,,^)2+  (fi^,2)] 

air 


1/2 


(3) 


dQair    ,,dm  2    d  Cpair  2         ,d <Tout-Tin>  2 


1/2 


air 


=  ((— )   +  ( 


)   +  (- 


m 


pai  r 


T    -T. 

out   in 


■)  ) 


-      dQ  .   2     dADr  2     ,._  2 

(4)   S-  ((^1I)   ^  (JT^)       +  (^)  > 
h      gair       APL 


1/2 


AT 


dNuv^  ^TT  o  dD,  ,  2     dK  .  _  2 

(5)     _!l£  =  ((^i)2  +  (]rM)  +  (irail)   ) 

Nu,  h  hd        air 
hd 


1/2 
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dRe 


(6) 


Re 


.  -       •  2     dD,  ,  2     du  .   2     dA   2 


1/2 


hd 


m 


hd 


ai  r 


The  following  quantities  had  their  uncertainty 
calculated  by  dividing  the  estimated  error  in  the  quantity 
by  the  value  of  the  quantity,  and  these  uncertainties  are 
considered  constant  for  the  range  of  values  in  this 
experiment : 


Quantity 
A 

Ac 

APL 

C 
pair 

Dc 


hd 


air 


atm 


id 


air 


Uncertainty 
.0093 
.0408 
.0080 
.0020 
.0400 
.0396 
.0030 
.0038 
.0010 
.0003 
.0080 
.0020 
.0054 


The  following  quantities  had  the  worst  case  error  for 
hydraulic  Reynolds  numbers  below  17000: 
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Quantity 

Tblk 

T. 
in 

orf 

T 
out 

T 
wo 

T    -T. 
out   in 

AT 


Error 


+  .1°C 


+  .1°C 


+.1°C 


+  .1°C 


+.1°C 


+  .2°C 


+  .2°C 


For  hydraulic  Reynolds  numbers  greater  than  17000  the 
error  in  the  outlet  temperature  was  much  higher,  as 
mentioned  previously,  and  the  following  errors  apply: 


Quanti  ty 
Tblk 


T. 

in 


orf 


out 


wo 


T   .  -T. 
out   in 

AT 


Error 


+  .2°C 


+  .1°C 


+  .1°C 


+  .3°C 


+  .1°C 


+  .4°C 


+  .3°C 


The  uncertainties  for  the  following  quantities  were 
calculated  using  the  sample  data  for  the  outer  curved  test 
section  at  a  hydraulic  Reynolds  number  of  15800: 
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Quantity 
Pl 

AP 


orf 

T   . -T. 

out   in 


AT 
P 


air 


m 


Q  • 
vair 

h 


Nu 


hd 


Re 


hd 


Uncertainty 
.0159 
.0108 
.0003 
.0521 
.0089 
.0159 
.0139 
.0540 
.0553 
.0681 
.0588 


The  uncertainty  for  the  following  quantities  were 
calculated  for  the  outer  curved  test  section  at  a  hydraulic 
Reynolds  number  of  21200: 


Quantity 
Pl 

AP 

orf 

T    -T. 

out   in 

AT 

p 
m 

Q  • 

vai  r 


air 


Uncertainty 
.0043 
.0053 
.0003 
.1039 
.0111 
.0043 
.0105 
.1044 
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h 

Nu 
Re 


hd 


hd 


.1053 
.1126 
.0581 


By  multiplying  the  quantity  by  its  uncertainty  the 
possible  error  can  be  calculated.   The  previous  uncertainty 
calculations  for  a  hydraulic  Reynolds  number  of  15800 
yielded  Nuhd  ■  40.07  +2.73  and  a  Rehd  =  15800  +  928.   For 


the  run  at  a  hydraulic  Reynolds  number  of  21200  the  Nu 

44.71  +  5.03  and  a  Re,  ,  =  21200  +  1232. 

hd         — 


hd 


The  major  source  of  uncertainty  in  the  average  Nusselt 
number  is  the  uncertainty  in  the  difference  between  the 
outlet  and  inlet  air  temperature.   If  this  uncertainty  can 
be  brought  below  four  percent,  then  the  uncertainty  of  the 
average  Nusselt  number  will  be  dictated  by  the  uncertainty 
of  the  measured  dimensions  of  the  channel.   In  particular, 
the  channel  height  (Dc)  which  has  an  uncertainty  dD C/D   = 
.0400.   The  uncertainty  of  the  channel  height  can  only  be 
reduced  by  increasing  the  channel  height,  and/or  using 
materials  with  better  tolerances. 
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